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New Format Presentations for Infrared Spectral Emittance Data 
I. INTRODUCTION 

Under the general title of "Supporting Research and Technology Effort 
in Identification of Geological Materials by Remote Infrared Spectroscopy" 
the Remote Sensing Laboratory of School of Earth Sciences, Stanford 
University, has been carrying out infrared (IR) radiance measurements from 
geological materials since 1967. These have involved laboratory and field 
spectroscopic measurements both on the ground and airborne. The net result 
of this work is a proven, feasible system for airborne use^over terrains 
with minimal vegetation. The perturbing effect of a long atmospheric path 
is still being evaluated and this may be a great barrier to successful 
spacecraft usage, although this in no way detracts from its highly successful 
operation from altitudes of several thousand feet. 

Recently parallel studies at University of Michigan £y Vincent and 

( 2 ) 

Thompson > based directly upon results of this work, but operating with 
image- forming airborne scanners, have shown that the spectral emittance con- 
cepts can be utilized in emittance-ratio imagery which broadly depicts the 
silicate composition of the terrain. This is a much more attractive 
technique than using spectral curves as the imaging mode of data presentation 
is much more meaningful to geological users. In addition the possibilities 
of analog processing the multichannel- imaged data are most attractive both in 
ease of data handling and cost. Considerable savings of time and funds can 
be made providing the requisite spectral resolution can be retained. A trade 
off exists between signal- to-noise (S/N) ratio and the spectral bandpass of 
each channel and channels which are too wide sacrifice spectral discrimination 
for high S/N ratios. 

- 1 - 

1. Lyon, R.J.P. (1972) "Infrared Spectral Emittance in Geological Mapping: 

Airborne Spectrometer Data from Pisgah .Crater, California"., Science , 175 , 983-986 

2. Vincent, R.K. , and Thompson, F.J. (1972) "Rock-type Discrimination from Ratioed 
Infrared Scanner Images of Pisgah Crater, California", Science , 175 , 986-988. 



II. OBJECTIVES 


Stanford University proposed three main objectives in the concluding 
phase of this study effort, 

A. Calculation of Emittance Ratios (Task 2.2) 

A first task involving the new concept of "emittance-ratio images" . 
These are important because they directly correlate with the surface 
composition of the terrain in an imaged form . They are obtained from the 
analog signals of a multichannel scanner which has several channels in the 
thermal infrared region. Such data could be readily obtained from the 
channels No. 17-21 of the MSDS scanner being prepared for the C-130 at MSC. 

By "banding together" the data in our higher-resolution airborne spectra, 
we have attempted to predict the response of the scanner over these terrains. 

B. Comparison of IR Spectral Emittance Data with K-Band Scatterometer 

Data (Task 2.3 ) 

The second task followed initial research efforts in comparing 
the airborne scatterometer data over Pisgah Crater, with the infrared spectra. 
Both methods produce "spectral curves" (backscatter versus look-angle, and 
emittance spectra) but represent very different skin depths for the information 
content in their data. 

The abstract (included in the appendix) is from a paper presented 
at the Seventh International Symposium in Remote Sensing of Environment. This 

2 

paper reported the use of a color-bar generation system (Digicol-I S) to 
generate a given color for each input IR spectrum. Prior to this stage of 
the analysis, the BMD07M stepwise discrimination program was used to identify 
the 3 wavelengths for best "separation" of the materials. These 3 emittance 
values were quantized within 16-level steps, and used to pre-set the 16-level 
matricies for the red, blue and green color guns of a color TV monitor. 
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The unit could handle 16 such matricies at once as vertical color bars across 


the TV tube. The eye then could rapidly identify similarities and contrasts 
in the data sets. Due to the familiar problems (color plate costs, etc.) in 
color reproduction the paper was not prepared in a final form for that 
Symposium volume. These techniques are a very interesting way to compare 
spectral (line-trace) data, and comparisons can be made with the IR and the 
X-band scatterometer data (5th Symposium volume by the University of Kansas 
group). Contrasts of skin depth such as appear with the windblown-sands 
covering the basalt flows, may be noted. 

C. Standard Infrared Spectral File (Task 2.4) 

The third task was to take all our spectra fro m both the ground and 
the airborne collection and prepare them in a standard format file. This 
is important for their use by any second party (who does not have our 
familiarity with their collection methods), and thus our accumulated expertise 
can now be passed on. These spectra have been collected from geologically- 
selected targets and 75-100 sets have been compiled. 

Where possible spectra represent rocks in their natural field 
location. In some sites (Sonora Pass) this was no longer possible and 
specimen rocks have been collected from sites. 
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III. RESULTS 

A. TASK 2.2 CALCULATION OF EMITTANCE- RATIOS 

1. Existing Spectral Data (MX108) - "Science” Article 

Over 4300 spectra were collected on magnetic tape from the 4 flight 
lines over Site 2 (Pisgah) on MX108. A total of 514 of these have been 
segregated into geologically-significant (and differing) categories and 
analyzed. Rather than repeat a description that is already in print, we have 
included the Science article verbatim. 


Infrared Spectral Emittance in Geological Mapping: 

Airborne Spectrometer Data from Pisgah Crater, California 

Abstract. Measurements of spectral omittance in the infrared region from 6.8 to 
13.3 micrometers were made with an airborne spectrometer at a rate of six spectra 
per second, on flights 650 meters above the olivine basalt flows at Pisgah Crater in 
the southern Californian desert. The spectra show chemical and mineralogicai 
diQerences that can be related to differences in the terrain below the aircraft . 


At a recent symposium (l, 2), two 
presentations were made of the tech- 
niques of independently performed geo- 
logical mapping from the air, over the 
same terrain. This report emphasizes 
the nonimaging technique involving 
spectral measurements in the infrared 
from 6.8 to 13.3 pm taken while flying 
at 650 m (2000 feet) above the olivine 
basalt flows and alluvium and dry lake 
beds (Lavic Lake) at Pisgah Crater, 
near Barstow, California. The instru- 
mentation (J), includes a boresight 
camera by which one can relocate the 
ground track of the (7 mrad, 0.4°) 
circular field of view of the spectrom- 
eter (4). 

A total of 514 of the 4300 spectra 
were collected (4) in four flight lines 
of total length 28 km (5). They were 
separated into 31 geological groupings, 
located as black and white bands on the 
large-scale photographs (Fig. 1) taken 

Reprinted from 

3 March 1972, Volume 175, pp. 983-986 


at the same time along the flight line 
F-F* (6). 

The raw radiance spectra from the 
rock and soil surfaces were ratioed by 
an average “water body” spectrum 
(average of 50 spectra) obtained by 
flying at the same altitude over a nearby 
lake (7). The emittance spectra were 
then inverted (5) and normalized for 
statistical studies, by setting their means 
at 0.0 and their standard deviation 
(S.D.) at 1.0. By this transform all the 
spectra have the same amplitude range, 
which permits more precise compari- 
son of their information content. Care 
should be taken in using these normal- 
ized spectra, as they are no longer 
numerically the same as absolute emit- 
tance (used for calculating tempera- 
tures from radiance levels). 

These selected spectra are then ana- 
lyzed in two formats. In the first, emit- 
tance values for 10 to 50 spectra [mean 
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Fig. 1. Photograph of Pisgah Crater (PC), center, and three image strips (10); 1 mile = 1.6 km. Fig. 2. Mean spectral emittance curves ( 11 ) 
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and ± 1 S.D.; (9)] are plotted as a 
function of wavelength. In the second, 
the standard deviations for 10 to 50 
spectra are plotted in the same way. 
The standard deviation plot provides a 
rapid estimate of the variability of the 
data. It generally has one of two shapes, 
low and flat if there is little variability 
in radiance, or higher and bulging if 
there is considerable variability. An in- 
teresting standard deviation plot for the 
dry lake sediments of Lavic. Lake 
(locality 28) shows a third type, with 
a pronounced upward bulge only where 
the reststrahlen effect {1-3, 10) indi- 
cated a chemical and mineralogical 
variation in an otherwise constant lake- 
floor terrain, which is clearly of geo- 
logical significance. 

The means of the spectral groups 
(or “mean spectra”) may be compared 
by visual inspection, for example, by 
direct overlay of tracings (Fig. 2), or 
group populations may be studied by 
discriminant analysis (77). Several sub- 
groups have been selected on the basis 
of more subtle features (such as weather- 
ing and surface chemical variability). 
In Fig. 2, B to H, the mean spectra of 
ail the subgroups arc plotted as single 
curves to show their similarities within 
a geological class. Thus, each spectrum 
represents the average of many individ- 
ual spectra sequentially observed over 
the distance indicated by the bars on 
F-F f in Fig. 1. 

The geology of the Recent lavas, 
rocks, and soils has been described 
{12, 75). The spectra can be correlated 
with geology and meaningful, system- 
atic variations appear in the flight data. 
The spectra of the “younger” alluvium 
(Fig. 2B) and “older” alluvium (Fig. 
2C) groups are similar, with a single, 
strong sloping minimum at 9.1 to 9.2 
fun. Detailed examination shows that 
the pattern in Fig. 2B is displaced to 
shorter wavelengths and that the shoul- 
der at 9.5 /on is absent, which indicates 
a higher quartz (sand) content in the 
younger materials (72). Comparable 
similarities arc shown by the three 
olivine basalt flow types (Fig. 2, F to 
H), all of which show a single, sharp 
minimum at 9.45 to 9.55 jum. The dif- 
ferences between the basalt spectral 
types were emphasized by separating 
the subgroups that showed a weak mini- 
mum at 10.97 /im (Fig. 2, F and G) 
from those that had a pronounced fea- 
ture there (Fig. 2H). A further separa- 
tion was made by means of the broad 
pattern around 11.5 to 12.0 /tm (in 
Fig. 2, F is flatter there than G). Com- 


parable spectra of (polished) granodi- 
orite (Fig. 2A) and gabbro (Fig. 21) 
specimens have been included to show 
the closer similarity of the basalt to 
the gabbo spectrum (74). 

The most interesting group (Fig. 2D) 
represents spectra from areas (3, 6, 
and 7 in Fig. 1) where blown sand 
now rests (patchily) in depths greater 
than the optical depth for these sili- 
cates. Thus, where cover is complete 
the spectrum of sand (here equivalent 
to the younger alluvium. Fig. 2B) 
should appear; where it is not com- 
plete the basalt spectrum should be 
evident. Within the group in Fig. 2D, 
the spectral mean IIC (locality 3) shows 
the younger alluvium pattern, while IIA 
(locality 7) is most like the basalt spec- 
tra (type 2, Fig. 2F), which establishes 
the ability of the airborne system to 
discern variations in rock composition. 
This variability is also clear on the 
A-B section of the ratio-imagery in Fig. 
i U0). 

The spectra of dry lake sediments 
(Fig. 2E) from Lavic Lake present an 
enigma. Although creamy white in 
color, these fine-grained clays (72) 
consistently yield spectra (75) similar 
to those of the type-2 basalt flows (see 
Fig. 2B and localities 27, 28, and 29 
in Fig. 1). This was true in 1965 in 
ground measurements along the same 
line. The new airborne data now sup- 
port the earlier (still unexplained) find- 
ings. Similar support is gained from the 
ratio-imagery of the same area observed 
by Vincent and Thomson {10). There- 
fore, three pieces of evidence point to 
the similarity of the clay spectra to 
those of the nearby olivine basalt. 

In summary, (i) infrared emittance 
spectra taken from the air over geo- 
logically selected areas across the Pisgah 
Crater lava flows show similarities with- 
in (and contrasts between) the areas. 
The spectral differences can be used to 
separate the flows, (ii) Within the lava 
flows themselves spectral types can be 
defined that (at the moment) are not 
clearly related to the mapped flow stages 
(flows I, II, and III); that is, the 
spectra are subtly depicting some other 
parameter than that used to differentiate 
the three flows in the field and on aerial 
photographs, (iii) Windblown sand on 
basalt shows spectra of sand, but where 
sand patchily covers the flow in one 
resolution cell the spectra of both the 
sand and the basalt appear in a com- 
posite pattern, (iv) The spectra of the 
dry lake sediments of argillic silts super- 
ficially resemble those of the basalt 


flows, but group variability in chemical 
and mineralogical composition is shown 
by the shapes of the standard deviation 
plots, (v) Each of these spectral simi- 
larities and differences may also be ob- 
served in the imagery prepared by ratio- 
ing concurrent radiance levels in two 
adjacent wavelength channels 2.5 /im 
wide (10). This is significant, as im- 
agery is more practical to use than 
spectral curves. What cannot be ex- 
plained yet is that this occurred with 
overlapping band-pass filters of such 
width (8.1 to 10.9 /Am and 9.4 to 12.0 
/un). These bands (Fig. 2) must rep- 
resent the integration of all the spectral 
information within their bounds and 
express it as an average value, rather 
than show all the finer points of spec- 
tral differences evident on the curves. 
A more precise separation of rock types 
can be effected by using nonoverlapping 
or narrower bands, or both, even with 
the lowered signal-to-noise ratio in- 
curred by the lessened energy through- 
put. Such a multichannel system as that 
being built for the National Aeronautics 
and Space Administration’s aircraft pro- 
gram (76) can become a geological 
mapping tool. 

R. J. P. Lyon 
School of Earth Sciences , 

Stanford University , 

Stanford , California 94305 
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follows, “play a surface (dry lake sediments) 
is a hard dense compact argillic crust con- 
sisting of approximately 79 % clay, 20% gran- 
ular components, 0.2% accessory minerals 
and a trace of saline minerals.” 

13. J. D. Friedman [C/.S. Geol. Surv . Tech. Lett. 


NASA-20 (1966), p. 4] gives the composi- 
tion of the basalt flows at Pisgah Crater: 

(flow II) total feldspar 15.4 percent and total 
ferromagnesian 9.4 percent; (flow HI) total 
feldspar 38.8 percent and total ferromag- 
nesian 5.1 percent. Flow HI is more 

spathic than flow II, at least at the two 

points sampled near the crater. Herein lies 
one of the major problems in relating re- 

gionally variable (airborne) composition data 
to those of classical geological studies, which 
are usually from selected points. 

14. See note 3 in (7). Standard rock I: gabbro, 
contains plagioclase (60 percent anorthitc 
molecule content), augite, and a little biotite; 
standard rock A, granodiorite, contains bio* 
tite, quartz, epidote, and plagioclase with 
orthoclase. 

15. The marked drop-off near 8 Am (In Fig. 2E) 
or the correspondingly high maxima at 8.8 
Am may be due to the Christiansen effect in 
these fine-grained materials. This, however, 
does not fit for the lavas (Fig. 2, F and G). 

16. The 24-channel scanner [E. M. Zaltzeff, C. L. 
Korb, C. L. Wilson, Inst. Electr. Electron . 
Eng. Trans. GE-9, 114 (1971)] has six chan- 
nels selected within the thermal band, as 
CH 16 (6.0 to 7.0 Am), CH 17 (8.3 to 8.8 
Am), CII 18 (8.8 to 9.3 Am), CH 19 (9.3 to 
9.8 Am), CH 20 (10.1 to 11.0 Am), CH 21 
(11.0 to 12.0 Am), and CH 22 (12.0 to 13.0 
Am). The data in (i) would represent the 
combination of channels 17 to 20, ratioed 
with the combined channels 19 to 21. My 
spectral data (Fig. 2) indicate that the Pisgah 
geology would be more clearly defined by 
using channels 17, 18, and 19 (either singly 
or combined), ratioed to channel 20. Channel 
21 would still show some effect of chemical 
compositions (particularly In femic rocks). 

17. This research was supported in its entirety 
by NASA contract NAS9-7313 with NASA/ 
MSC, Houston. This financial support is grate- 
fully acknowledged. This is Technical Report 
No. 71-5, Remote Sensing Laboratory, Stan- 
ford University. 


21 September 1971, revised 22 November 1971 


2. State of Art at End of 1971 Contract 

The state of processing represented by the above article has been 
materially advanced during this new period, but the details will be covered' 
in the next section. 

Again, the material is best covered by reproducing a published paper 
by Lyon and Marshall, 1971, in IEEE Transactions on Geoscience Electronics. 
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Operational Calibration of an Airborne Infrared 
Spectrometer Over Geologically 
Significant T err ains 

R. J. P. LYON and A. A. MARSHALL 


Abstract — A three-instrument infrared spectral emittance experi- 
ment, comprising a rapid-scan spectrometer (6.7-13.3 ^m), radiome- 
ter (10.375-12.1 mto), and boresight camera, has been flight tested 
over selected geological terrain in central and southern California 
and Nevada. Pre- and post-flight calibrations of the infrared spec- 
trometer were performed both by using polished samples of “stan- 
dard"' rocks (quartz diorite and gabbro) as well as the more familiar 
blackbody radiance standards. From these latter spectra the instru- 
ment transfer function (/W\) was derived. In-flight calibrations of 
wavelength were achieved by the rapid insertion and removal of a 
polystyrene film in the optical train of the spectrometer, as poly- 
styrene is a material whose transmission spectrum is constant and 
well known. By flying over a body of water (<*. 7 -u. 3 *m = 0.98) and 
recording the radiance spectrum of that target one can determine the 
transmission spectrum of the atmospheric path between the aircraft 
and the water (at least to a first approximation) as both the spectral 
emittance of lake water and the optical transfer functions of the in- 
strument are known or can be calculated. So far, flights have been 
made only at low altitudes (2000 ft above the lake), with the lake 
surface at 2000 ft (near Pisgah Crater, S. Calif.) or 6000 ft (Mono 
Lake, E. Calif.) above sea level. The lake should be in the area to be 
studied geologically. If the flight altitudes over the study areas are 
consistent with those over the lake, then the effect of the airpath can 
be evaluated relative to the spectral information from the geological 
targets. 

Introduction 

T O ACCURATELY deduce the surface tempera- 
ture of terrain from its infrared brightness tem- 
perature, the surface emittance as well as the 
background radiant emittance must be kn<?wn. Gen- 
erally when such measurements are made around S-14 
/xm a simplifying assumption is used that 1.00, 

i.e., the target is a blackbody. Sometimes a graybody 
emittance of 0.9 or similar factor is assumed to ease the 
conscience of the researcher involved. One method of 
arriving at the emittance integrated across a given 
passbaiul is to measure the spectral emittance and then 
integrate for an averaged value. This experiment de- 
scribes equipment and analytical techniques by which 
this may be achieved. 

A geologically significant more sophisticated experi- 
ment which relies upon the non black body behavior of 
silicate rock materials typically making up planetary 
terrain, is to use the spectral emittance in this band to 
derive the chemical composition of the terrain being 

Manuscript received February 26, 1071; revised March 24, 1971. 
This work was supported by NASA, under Contract NASO-7313. 

The authors are with the Remote Sensing Laboratory, School of 
Earth Science, Stanford University, Stanford, Calif. 04305. 


overflown [l ]— [3 ]. Subsequently, one can integrate the 
spectral data to obtain a suitable wide-band value. The 
background information on this method is detailed in 
several references, and the method has been reduced to 
practice both in the field and in airborne measurements 

[31. 

In the airborne mode, infrared radiance spectra (with 
a resolution X/AX = 1O0) were taken six times a second 
over the bandpass 6.7 to 13.3 /xm. With every second 
spectrum a 35-mm boresight photograph was taken to 
locate the precise ground-track of the sensor system 
after the data flight. Table I lists the equipment charac- 
teristics. 

Also included as the third instrument in the infrared 
“pallet” was an infrared radiometer, a relatively broad- 
band sensor (10.375 to 12.1 jxm), which has its bandpass 
centered in a region of high atmospheric transmittance. 
This band exhibits consistently high terrain emittance 
and does not show the spectral departures from a black- 
body on which the geological experiment is based [3]. 
The radiometer served as a monitor to ensure that there 
were no marked temperature changes over the target 
which might be mistaken for spectral emittance 
changes. Within a given spectrum, temperature changes 
of less than 1°C were allowed — if higher than that level 
the spectrum was rejected. 

Calibration Concepts 

The operational calibration was of three main types, 
pre-flight, inflight, and post-flight; the pre- and post- 
flight sets being performed immediately prior to take-off 
and after landing. 1 The aircraft engines usually were not 
running and the internal systems were connected to 
ground power (sometimes the auxiliary power unit 
(APU) was used). There could be differences between 
the two ground power sources (and with the aircraft 
engine sources themselves) but the most significant 
“noise” problem at this stage is the degree of mechanical 
vibration generated by the APU, and the engines which 
induces microphonic signals into the infrared detectors. 

The in-flight calibrations are more simple, being 
restricted to a wavelength check performed by inserting 

1 Cryogenic hold-times were in excess of eiijht hours; thus long 
flights could be made and post-flight calibrations obtained without 
refilling the cryogenic supplies. 
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TABLE 1 
Infrared Pallet 


Airborne Rapid Scan Spectro meter 

Scan wavelength 

6.76-15.30 /m i with 100 elements per S{)ec- 
trutn. The CVF* wheel has 2 similar spectral 
octaves — one from 0° to 180°, and one from 
180° to 360° 

Scan period 

0.150 s (6 spectra /s) 

Field of view 

0.4 degree square (7 mrad) 

Detector 

Hg-dopcd germanium, time constant less 
than 1 /is, cooled by liquid helium 

Essential output 
• signals (four) 

a) spectral radiance output (analog) 

b) wavelength ramp (analog, not presently 
used) 

c) wavelength (peripheral-edge coding) 
pulses, every 2°, or 90 j>cr S|>ectntni, 180 
per rotation of the CVF (See Table 11) 

d) a spike pulse, (at 0°) was used to fire the 
bo resight camera (used for location pur- 
poses) 

Accuracy required 

10-bii, i.e., better than 0.1 percent 

Infrared Radiometer 


Filter bandpass; 
sampling frequency 

10.375 to 12.1 «m approximately 60 tempera- 
ture measurements [>er second, i.e., ten to 
every s(>ectrum (or 1 every 9 spectral ele- 
ments) 

Field of view 

0.4 degree, circular (7 mrad) 

Detector 

Hg-doped germanium, time constant less 
than 1 /is, liquid helium cooled 

Essential output 
signals (one) 

radiance signal sampled 60 times/s (ana- 
log) 

Accuracy required 

10 bit, i.e., better than 0.1 percent 

Boresight Camera 


Type 

35-nun framing camera, with film-recorded 
clock and frame counter, electrically pulsed 
by output command from sj>ectrometcr (at 
approximate rate 3 s) 

Field of view 

approximately 5° to yield telescopic view of 
the target. Camera pulse originates 5 ms be- 
fore the no. 1 data pulse, i.e., just past the 
0° jjosition 


• CVF: circular variable filter, a circular dispersive element. 


a polystyrene film into the spectrometer optical train. 
The known spectral transmission of polystyrene pro- 
vides a characteristic signal readily observed on the 
data tapes. It was also used as a signal for •Mine-start” 
and “line-stop” in the later data analysis steps. Fig. 1 
shows a group (A* = 5) of airborne polystyrene spectra. 
From this, Table 11 was verified. 

A more fundamental type of in-flight calibration was 
performed by th ing the infrared system over a body of 
water. Water has a well-known spectral omittance vary- 
ing only slightly from a grayhotly value of 0.98 [6 j over 
this bandpass. Thus a radiometric check between 
spectrometer and radiometer may be made in the data 
analysis steps and brightness temperatures compared. 
As will be seen later, the AIKIWTII spectral transmit- 
tance between the water and the aircraft may be defined 
by this aspect of the flight program. 



Fig. 1. Average and ±<r of AT = 5 polystyrene transmission spectra 
airborne over Pisgah target. A f = 5. Times from 17 : 13 : 14.442- 
17: 13: 16.000 GMT (local time PST = GMT-g h). Handwritten 
figures refer to the precisely known absorption lines for poly- 
styrene. 


TABLE II 


Counter Point (CP) Versus Wavelength* 


CF 

X (/on) 

CP 

X (pm) 

CP 

X ( jim ) 

91 

6.760 

121 

8.965 

151 

11 .380 

92 

6.825 

122 

9.045 

152 

11.445 

93 

6.890 

123 

9.135 

153 

1 1 . 525 

94 

6.960 

124 

9.205 

154 

11.600 

95 

7.035 

125 

9.270 

155 

11.680 

96 

7.085 

126 

9.370 

156 

11.750 

. 97 

7.160 

127 

9.455 

157 

11.835 

98 

7.225 

128 

9.530 

158 

11.895 

99 

7.290 

129 

9,615 

159 

11.970 

100 

7.365 

130 

9.690 

160 

12.040 

101 

7.4+5 

131 

9.775 

161 

12.120 

102 

7.510 

132 

9.865 

162 

12.190 

103 

7.585 

133 

9.960 

163 

12.260 

104 

7.660 

134 

10.010 

164 

12.330 

105 

7.750 

1 35 

10.100 

165 

12.405 

106 

7.810 

136 

10.170 

166 

12.470 

107 

7.880 

137 

10.255 

167 

12.540 

108 

7.955 

138 

10.340 

168 

12.605 

109 

8.025 

139 

10.420 

169 

12.670 

110 

8.105 

140 

10.490 

170 

12.745 

111 

8.185 

141 

10.580 

171 

12.815 

112 

8.265 

142 

10.660 

172 

12.895 

113 

8.325 

143 

10.740 

173 

12.955 

114 

8.410 

344 

10.810 

174 

13.020 

115 

8.485 

145 

10.900 

175 

13.085 

116 

8.565 

146 

10.970 

176 

13.150 

117 

8.640 

147 

11.070 

177 

. 13.220 

118 

8.720 

148 

11.150 

178 

13.285 

119 

8.800 

149 

11.220 



120 

8.885 

150 

11 .305 




s Original data sheet for points 0-90 and 9 1 -ISO for the 1908-1970 
ARSS CVF wheel, provided by NASA MSC by letter TF2 LE 120- 
68, March 25, 1968. To make the data covered by both halves of the 
wheel more equal 1 CP is dropped from the leading data set for each 
>itlc. Thus M.SC CP 92 becomes Stanford CP 91, 179 becomes 178, 
etc. Operational verification of these values may be made by the 
sjiectrum of polystyrene in Fig. 1. 


Pre- and Post-Flight Calibrations 

In addition to the above, there were normal engineer- 
ing-type calibrations of tape recorder channels, etc., of 
the type usually performed iu all high quality data 
gathering operations [4], The calibrations described 
here refer to geological and meteorological calibrations 
which are designed to recover known spectral charac- 



LYON AND MARSHALL: CALIBRATION OF AIRBORNE INFRARED SPECTROMETER 


10 


teristics from the data tapes as an overall check of the 
system. In detail these arc the following. 

1) Observation of an external blackbcdy source lying 
on the runway beneath the stationary aircraft. This 
blackbody was temperature controlled at 40°C (313°K), 
and is called 40°C BB-EXT.* 

2) Sequential observation of a pair of “standardized” 
rock specimens (rock A, and rock B), lying on the run- 
way beneath the stationary aircraft. Those were heated 
in an oven to approximately 40°C (313°K). The speci- 
mens were 20 by 20 cm across and 2.5 cm thick, thus 
having considerable thermal inertia. 3 Rock A was a dark 
gray to medium gray gabbro, and rock B was a light 
silver-gray granodiorite, two chemically, and mineral- 
ogically, distinct igneous rocks. The front surfaces of 
these “standards” were highly polished to decrease 
scattering effects from their front surface making it a 
more precise spectral emissivitv standard. 

In order to obtain high contrast in the spectral signa- 
ture data, i.e., a large signal differential for small spec- 
tral emissivitv variations, the background radiation 
reflected by the sample into the instrument must be 
small compared to the emitted radiation. This criterion 
was achieved by positioning the highly polished sample 
sufficiently below the aircraft and angled appropriately 
to ensure that “cold” sky was the effective background. 4 
Since all pre- and post- flight calibration measurements 
were performed under moderately low humidity and 
cloudless conditions, the sky radiation (between S jum 


* Extended Source Rlackhtxly, Barnes Engineering, Inc., tempera- 
ture controlled to ±0.5 J C. 

1 The rocks were selected by a geologist from tombstone material 
(because of the overall high polish readily available across large slab 
areas). Thin sections were made and the correct rock type defined by 
optical means. The slabs were cut into two, and one half of each 
retained by Stanford University Infrared Lab. to calibrate our other 
spectrometer and data recording systems. Rock A: Gabbro, contains 
plagioclase (An so), augitc. and little biotite, Rock B: Granodiorite, 
contains bictite, quartz, epidote, and piagioclase, with orthoclasc. 

4 The spectrometer measures IIxmo. Assume a target which en- 
tirely fills the instrument field of view. 

= H"xa(l — r\«) + rxofSuIUiX 4- — Six)] 

where 

Uxtefo effective radiant emittance seen by the radiometer at 
wavelength X; 

W\ a radiant emittance of the air at wavelength X; 

T\ a atmospheric transmission at X; 

Six target emissivitv at X; 

H'rx radiant emittance of the target at X; 

H\b irradiance from I he background (sky, terrain, mom walls, 
or ceiling, etc.) incident on the target at wavelength X. 

Note t he following: 

1) H x «. r\ 0 , ir tX , and U\u are integrated over the s|>ectral 
bandpass of the instrument. The difference between the integrated 
values and monochromatic values at wavelength \ depends on the 
siicctral bundwidths of the atmospheric absorption /emission band?-, 
the bandwidth of the spectral features on the target, and similar 
considerations for the background relative to the instantaneous 
sjvctral bandwidth of the instrument. In the 7-M-jim region these 
differences are generally small for AX/A <0.02. 

2) /A* depends upon r x „, 1! Xa , sky conditions, terrain temperature, 
the target angular asjiect, and instrument view angle relative to the 
zenith. If U’ Xr = //*« no sjcciral signature of the target is discernible. 
IIxb can be better controlled or minimized with a polished (specular) 
sample. 

3) rxo has practical limits between zero and one and has both fine 
and coarse sj>cciral features. 



Fig. 2. Calculated radiance difference (W-cm - * sr -1 Mm” 1 ) for a 
blackbody exterior target at 40°C (313°K) and an internal black- 
body reference at 60 3 C (333 °K). Wavelength counter points 
appear as the abscissa, from 91 to 178 (from 6.8 to 13.3 Mni), and 
are common to all graphs. 


and 13.5 fim) was infinitesimal compared to that of the 
40°C rock specimen and was neglected from further 
consideration. Identical procedures were followed dur- 
ing pre- and post-flight calibrations. 


Transfer Function 04 t x ) Calculation 

The optical spectral transfer function (r x ) for the 
spectrometer may be derived jointly with the system 
gain A as a product form, which we simply call the 
spectral “transfer function” (zlr x ). The data processing 
steps outlined in what follows agree closely with those 
published for the Purdue (LARS) procedures originally 
performed for field calibration of BIock-Michelson inter- 
ferometer spectra [5]. 

Calculate for each wavelength step the radiance 
spectrum (Planck Law) for the 60°C internal blackbody 
reference (BB-INT). Calculate the radiance spectrum 
for the temperature controlled external blackbody 
(BB-EXT). Derive the radiance difference spectrum 
(TRUE). 

Fig. 2 shows this calculated difference 

iVXxRUE = -YAbb-int — ^Abb-ext (1) 


where iVX = radiance at wavelength A (jxm ). 

This difference A t X T rce can be related to the actual 
signal observed from the spectrometer (proportional to 
the difference between the two blackbodies, but modified 
by the spectral transfer function). Figs. 3 and 4 show the 
actual signal (OBSERVED) from an external black- 
body at 40°C (313°K); an average of 22 spectra with 
± lor limits indicated. The standard deviation levels arc 
a measure of the “noise” riding on the signal. From these 
data (Figs. 3 and 4) the spectral translcr function (.4r x ) 
was calculated (Fig. 5) 


SPECTRAN (.4r x ) = 


OBSERVED 


( 2 ) 


TRUE 
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Fig. 3. Observed radiance dilTerence (millivolts) for the same 40°C 
target and 60'C reference. A ' = 22 spectra for the average calcula- 
tion, ±l<r limits shown, X -668 mV. Noise (mean sigma, a) 
=s 13.5 mV. Pre-flight data, mission 108, day 1, times from 
15:17: 53.61 7-15:1 8 :00. 759 GMT. 



Fig. 5. Optical transfer function (.-In), from pre-flight 40°C BB- 
EXT data. Ratio calculated from Figs. 2 and 3. 

The inverse of the transfer function 1 /At\ or IN- 
STRAN was more useful in the calculations and this 
function was stored for data processing (Fig. 6). The 
principal straight-line component of the gradient is 
related to the wavelength sensitivity of the Ge:IIg de- 
tector, the system gain (.-1), and the optical trans- 
mittance (tx). Higher frequencies in the function are 
caused by the variations in tv exhibited by the CAT. 
Additional small effects of microphonia noise and inter- 
ference may account for some of the observed high fre- 
quency variations. 

AIRPATH Calculation 

When airborne spectra are collected while flying over 
a water body (Figs. 7 and S) a mean spectrum can be 
calculated which closely resembles the Maekbody differ- 
ence spectrum, OBSERVED (Fig. 3). Absolute levels 
(in volts) arc higher, as most water bodies are cooler 
than the 40°C external blackbody used in the pre-flight 
ground calibrations. Standard deviation levels of 17-20 
mV rins (Table 111) again represent noise levels but are 
not significantly higher than on the ground, a marked 
change from earlier flights with this system which 



Fig. 4. Comparable post-flight data after return. iV-25, a -19.1 
mV, X-692 mV (re., temperatures were not exactly at 40 C C). 
Times from 21 :21 :23.167-21 :21 :30.766 GMT. 



showed large microphonic noise. (Missions 56, 78, flown 
in 1967 and 1968. ) 3 

The atmospheric transmittance was determined over 
the two lakes observed in this flight (Palmdale Lake, 96 
mi west, and Shallow Lake, 25 mi north of the Pisgah 
lava field test site, which is at latitude 34.7°N, longitude 
116.4°W). They are a considerable distance apart, and 
somewhat distant from the test site, but it is difficult to 
find water bodies in the southern Californian deserts! 
One is forced, therefore, to rely upon the assumption 
that the airmass over the test site is the same as that 
over the lakes. In most areas where water bodies are 
more common and hence closer to the site this is a quite 
reasonable assumption. Here one may argue otherwise, 
but no simple operational alternative exists. 

AIRPATH is calculated from the airborne data over 
the lake by the following steps. 

1) Calculate the mean and ± ler spectra for the 
(LAKE) airborne data. Calculate lake brightness tem- 

5 OjHirational rins voltage measurements (“noise”) were made on 
the l R sfjectrometer output (with theCVF wheel stopped), as follows. 
MX l OS, flight t. Pre-flight — on ground 20 mV. 1915 GMT- — over 
water 25 mV. MX 108, flight 3. Over Mono Lake, 1641 GMT 24 mV. 
Over Mono Lake, 1940 GMT 24 mV. 
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52.400 GMT. mi— ion 108. flight I. Ocwpoint— 6 C C: air tempera- 
ture 23°C; RH 14 percent; radiometer Tl{R) =30.5°C. 


\ 



Fig._ 8. Observed radiance spectrum from Shallow Lake, iV = 26, 
ff = 17.7 mV, = 1000 mV. Brightness temperature of lake 
2^=3 0 C C. Times from 19:16:20/898-19:16:28.470 GMT. mis- 
sion 108, flight 1. Dewpoint— 6°C; air temperature 27°C; RH 
10 percent; radiometer Tl(R) = 33 °C. 1 


TABLE III 

Tabulated Reduced Data 


GMT Time 

Stanford/MSC Mission-108, Flight l (Pisgah, Site 2), Oct. 8, 1969 

N Signal Mean Std. Deviation* 

(mV) (mV) 

Rel. Error 6 

Temperature by 
Spectrometer 


Pre- Flight 






1518 

40°C BILE XT 

22 

669 

13 

0.0043 

40*0 

1516 

rock A 

30 

569 

15 

0.0048 

43 

1517 

rock B 

30 

595 

15(20)' 

0.0046 

42 


In-Flight 






1712 

Palmdale lake 

30 

1010 

15 

0.0027 

30°C 

1915 

Shallow lake 

30 

1000 

18(25)' 

0.0033 

30 


l’ost-Pliuht 






2121 

40'C BB-EXT 

25 

693 

19 

0.0055 

39°C 

2120 

rock A 

36 

512 

22 

0.0072 

45 

2120 

rock B 

28 

568 

24 

0.0080 

43 


Stanford/MSC Mission-108, Flight 3 (Mono, Site 3), Oct 10, 1969 





N 

Signal Mean 

Std. Deviation 

Rel. Error 

Temperature by 

GMT Time 



(mV) 

(mV) 


Spectrometer 


Pre-flight 






1423 

40°C BB-EXT 

30 

664 

16 

0.0042 

40°C° 

1422 

rock A 

30 

569 

15 

0.0048 

43 

1424 

rock B 

30 

593 

15 

0.0046 

42 


In-Flight 






1642 

Mono Lake run 401 

30 

987 

19(24)' 

0.0034 

30°C d 

1943 

Mono Lake run 402 

30* 

1080 

17(24)' 

0.0029 

27* 


• Standard Deviation (a) is defined as 


/(* - X) 1 

* " V x - i 

c By definition, and used as a calibration point for temperature. 
d West side. 


b Relative error is defined as 

E m YwW ’ 

• East side. 

1 Actual measured rms values in parentheses. See footnote 5. 


perature (7T), using the radiometer data, i.c., with 
broader bandpass. This unit was also recalibrated during 
the pre-ilight period [4]. 

2) Correct the mean lake spectrum for the known 
spectral cniittancc for water (e x ) [6]. 

3) Divide the correct (observed) lake spectrum by 
the product of a calculated blackbodv radiance spectrum 
for that temperature ( T /.) and SPECTRAX (-lr x ). 
This modifies the blackbodv radiance as though it had 
been observed through the instrument. If the lake had 
been at 40°C one could use the observed spectrum from 
the external 40°C blackbody. Generally, however, the 


lake temperatures are much lower (27-30°C), and hence 
their radiance must be calculated. 


expected radiance 

AIRPATH = — 

actual radiance 


LAKE calculated 
LAKE observed 


then 


AIRPATH 


(g box, at T° L ) X A r\ X €\ 
LAKEx 


(3) 


AIRPATH is an absorbance-like ratio which ranges 
from 0.84 to 0.96 with 10 to 12 sharp maxima. These are 







136 




generally in the positions of the atmospheric absorp- 
tion/emission bands, but may be more exactly caused 
by the high frequency components in the pre-flight cali- 
bration spectra which were ratioed to the lake spectra 
(Figs. 9-12). As might be expected, for a 2000-ft path 
in a low 6 humidity (desert) area, these absorptions are 
not strong between 7.5 and 13 pm. It is significant that 
both lakes surveyed this day showed very similar pat- 
terns with maxima at the same data points. On a sub- 
sequent flight (day 3) over Mono Lake several hundred 
miles to the north, both patterns for AIRPATH were 
again similar (Figs. 11 and 12), but differences exist be- 
tween the flight 1 and flight 3 pairs. 

Rock Standards — Operational Format 

The pre-flight spectra from the two polished rock 
standards are shown in Figs. 13 and 14 with the mean 
and l<r limits plotted. Immediate inspection reveals that 
they have different shapes as expected; rock A having 
two maxima at longer wavelengths than the more pro- 
nounced single peak for rock B. 

At first glance the curves appear to be the inverse of 
what one would expect from emittance data. A mo- 

1 Dewpoint— 6°C, aircraft air temperature +23°C. 


13 

IEEE TRANSACTIONS ON GEOSCIENCE ELECTRONICS, JULY 1971 X 



Fig. 10. AIRPATH, from Shallow Lake. Mean airpath (4p) =0.96, 
£**=* 0 . 017 . 



Fig. 12. AIRPATH, from Mono Lake 402 (?ite 3) r flight 3. mean 
airpath (Ap) =0.950, Oa p =0.013. A’ =30, 7 t=27°C, times from 
19:43 : 55. J 1 2—19:44 : 14.813 GMT, first 30 spectra used. Dew- 
point— 19°C; air temperature 10°C; RH 12 percent. 

ment's reflection, however, will show that if c x is lower at 
a given point it will appear to have a lower (colder) 
brightness temperature. As we are using a hot internal 
reference blackbody (60°C) these regions will show a 
greater radiance difference, and hence a higher millivolt 
level, for low c x regions. In effect, the raw plots are re- 
flectance not emittance spectra. 

For rapid chocking of the data either this form of 
presentation or its exact inverse (FLIPPED) will 
suffice. For most operational use we do not go through 
the calculations for absolute emittance detailed below, 
but merely “Hip” the data sets. All airborne spectra of 
the terrain are presented in tins inverted form for sub- 
sequent statistical processing [7 ]. 

Rock Standards— Absolute Emittance 

The absolute emittance from rocks A and B may be 
calculated by the following steps. 

1) Prepare a mean and ±1 cr average rock A spec- 
trum from the pre- and/or post-flight data (usually 
A 7 = 30). Determine the average brightness temperature 
from these spectral levels, e.g., X = 593 mV. Therefore, 
target brightness temperature was T A =42°C. 

2) Using the spectral transfer function (At\), trans- 
form the mean spectrum, wavelength interval by wave- 
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Fig. 13. Rock A mean sjicctnini for the gubbro slab at T =43°C, 
*Y = 56, from raw (i.e., in ^mouthed) pre-fli'Ja data, flight 1. Times 
are from 15: 16:42.190-15: 17:00.963 GMT. Teaks are seen at 
122 (626 mV) and 1.10 (660 mV) data points. Radiance levels at 
data point IDS would approximate that from a 4()°C blaekbody 
target, at MO would be that of a 25 C blaekbody. This is indica- 
tive of the changing sj)eciral cmitt nice. 




Fig. 14. Rock B mean spectrum for the granodiorite slab at 
r = 42°C, A r = 4 3, from raw pre-flight data, flight 1. Times are 
from 15: 17: 10.651-15: 17:23.674 GMT. Peaks are now at 114 
(643 mV) and 124 (673 mV), much shorter wavelengths than for 
rock A. 



Fig. 15. Radiance mean spectrum (±U) for rock A, pre-flight 
data, flight 1, for the above times (Fig. M): Ordinate is now in 
\Y cm -2 sr~ l jum“ l , after conversion from millivolts using opti- 
cal transfer function (.Irx). 7 a=43 : C, spectrometer. 

length interval, into a radiance spectrum, i.e., transform 
millivolts to \Y*cm~ 2 sr“ l /inr 1 . In a similar manner one 
can transform the standard deviation spectra as well. 
This is plotted in Figs. 15 and 16 for rock A and B as 
radiance spectra. 

3) Using the observed target brightness temperature 
(as above, of 42°C) calculate from the Planck Law. 
again one wavelength interval a t a time, the expected 
radiance for the target, assuming now that it is a black- 
body. 

4) Divide the “observed” radiance spectrum by the 
calculated blaekbody radiance spectrum to generate an 
emittance spectrum, Figs. 17 ancl 18. By the initial 
assumption of A' x average brightness temperature this 
yields emitlance values above 1.0. These then can be 
normalized to unity and the average “absolute” emit- 
tancc recalculated. 

Two items should be noted here. Inspection of Figs. 
15 and 16 reveals the very small departures from bl.ick- 
bodv behavior on which this geological experiment is 
based. Figs. 17 and IS are therefore magnifications of 
the fx ratios, and in this light the ± l<r variations in the 



Fig. 16. Radiance mean spectrum ( ± l<r) for rock B, pre- flight data, 
flight 1, for the above times (Fig. 14). Ordinate is now in W -cm -2 
sr~ l /im -1 . ra = 42°C, spectrometer. 

airborne data are strikingly small. As geologists rather 
than physicists we are not as concerned with the ab- 
solute levels in our data as with the homogeneity of the 
spectral data and its variability from one rock type to 
another. The strictly geological problems with this ex- 
periment are many and real, ancl added refinements in 
computational technique are premature at this stage. 

CONCI-C'SIOXS 

1) Optical transfer functions for a spectrometer have 
been calculated from operational field-type measure- 
ments. 

2) Standardized rock slabs with known spectral 
emittancc below 1.0 have been used in the pre- and post- 
flight calibrational steps in a routine manner. From 
these data one may ascertain if the spectrometer and 
data system arc functioning correctly, i.e., that a known 
spectrum can be recovered from the taped data record. 

5) In-flight calibrations of wavelength can be readily 
obtained by inserting polystyrene into the optical train. 

4) A tmusplierie . “absorption spectra” can be ob- 
tained in f he airborne mode by flying over lakes within 
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Fig. 17. Calculated “emit lance* ratio (or rock A target relative to 
a blackbody at (lie sann* average brightness tcni|)crature. Mean 
and ± l<r limits shown. Average omittance = t .00 by assumption, 
«ra« = 1.04 (108), «mi»=0.98 (144). When normalized so that the 
(max ^ l *00, then e n .m=0.94. 

or nearby the geological test site. These spectra have 
many of the characteristics of true atmospheric spectra 
of a short airpath (around 2000 ft) in a low humidity 
area. Additional work is necessary to confirm t his inter- 
pretation or to support the possibility that they are arti- 
facts of the ratioing of airborne data (which appear 
somewhat smoothed; to those of ground-based data 
which are more spiked. 

5) ln-flight spectral data over geological targets are 
being analyzed, and will be reported elsewhere [7 ]. Con- 
siderable careful study needs to be made, both of the 
airborne data as well as spectral omittance variabilities 
within a single rock type in the field. This should be 
tackled by ground -based mobile instrumentation as 
well as airborne to define the inhomogeneity of the basic 
target and separate this effect from artifacts of airborne 
measurements. Only then can we assess the “noise” un- 
certainties we presently perceive. 
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Fig. 18. Calculated “emittancc” ratio for rock B. Mean and ±1 <j 
limits shown. Average emit tancc = 1.00 by assumption, 

= 1.02 (106), (,,,10=0.97 (118). When normalized so that (mu 
= 1.00, then «„» ft =0.95. 
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3. Generation of Emittance Spectra - 1972 Procedures 

(a) Airborne Spectra 

The basic data obtained from the digital tapes is a set of voltages 
digitized at a corresponding set of wavelengths provided by NASA in previous 
calibration work. These voltages must then be converted with the aid of several 
calibration and normalization procedures into a form which resembles true 
emittance spectra. 

The voltage generated by the spectrometer, as a function of the 
spectral radiance incident on the spectrometer (N S p) > is dependent upon the 
temperature and thus the emitted spectral radiance of the internal reference 
blackbody N bb (T£ nt ) and the instrument transfer function At. 

V - At (N bb (T int ) - N gp )* (1) 

In general N S p is made up of two contributions 

(a) The radiance emitted by the surface as it is attenuated 
by the atmosphere (N sur x) and 

(b) The emittance of the atmosphere ((1-x) NjjhCTsky) assuming 
uniform composition and temperature distribution over the 
observing path length) . Thus the output volts given by 

V - At (N bb (T int ) - (1-x) N bb (T air ) - TN sur ) (2) 

To use the above equation three pieces of information are required to 
solve for N su r in terms of the output voltage. However, in studies involving 
relatively short paths (less than or equal to 2000 feet) and involving warm 
rocks (approximately 40° c) we have 

Ng U j- > N bb (Tatin) end (1— x) « x 
Thus we can simplify Equation 2 to give 

V = At (N bb (T int ) - X N sur ) (3) 

In the ground based calibration work at very short path lengths At and Xg were 
determined by the observation of two blackbodies , one at 40°C and the other at 
ambient temperature. Equation 3 could then be solved for the two functions At and xg. 


*A11 quantities mentioned are functions of wavelength. 



17 . 


Of the two 40° blackbody measurements made on the day of the 
measurements over Pisgah Crater, it was found that the pre-flight 40° black- 
body spectrum had a large amount of periodic noise (approximately 60 HZ) superimposed 
on the blackbody curve (see Fig.Al ). Until this was realized flight spectra 
calculated using this measurement as a calibration were found to be seriously 
perturbed, the effect of which was masked in the normalization process used 
at the time. Thus in all successive calculations the post-flight 40° black- 
body data were used in the instrument transfer function calibrations. 

With the knowledge of At it was then possible to calculate the spectral 
radiance reaching the spectrometer during the flights over Pisgah Crater and 
during the calibration flight over Shallow Lake. The spectrum of spectral 
radiances recorded off Shallow Lake were then examined and a radiant temperature 
computed for each of the wavelengths at which measurements were taken. The 
wavelength at which a maximum blackbody temperature (T^av) was found was then 
selected and the temperature of the Lake was then taken to be this radiant 
temperature. This is equivalent to assuming that at one of the wavelengths 
measured the atmosphere had a transmission equal to 1.0. The atmospheric 
transmission at all other wavelengths was then computed assuming that the 
radiant energy leaving the lake's surface at this wavelength was that of a 
blackbody of temperature I mav . 

These atmospheric transmissions were then used to compute the surface 
spectral radiance over Pisgah Crater. The wavelength dependence of this radiance 
was then examined and the wavelength within each recorded spectrum showing the 
maximum radiant temperature was again selected. The resulting radiant temperature 
was then assumed to be the rock temperature and the emissivity of the rock at 
the various wavelengths was formed by taking the ratio of the calculated surface 
radiance to the blackbody radiance at this maximum temperature. The previously 



Figure A1 

Preflight blackbody emittance spectra. 
Here the preflight calibration spectra 
have been treated in the same way as 
the rock samples. The upper curve is 
the Preflight 40°C Blackbody, the 
lower the Preflight Ambient Blackbody. 
The hot blackbody shows definite 
periodic noise which makes it unsuit- 
able for calibration purposes. 




selected groups of spectra were then separated and overall emittance spectrum 
for each geologic unit was computed, as well as the standard deviations at 
each wavelength, 

b. Ground Measured Spectra 

All the ground based measurements of rock emission spectra have been 
made with an Exotech Model 10 Spectrometer. The design of this instrument is 
essentially the same as that of spectrometer in the infrared pallet except that 
the reference blackbody is at ambient temperature in contrast to the 60° C black- 
body in the infrared pallet. 

Calibration measurements to determine the instrument transfer function 
were made using a water heated' blackbody source while wavelength calibration 
was obtained with polystyrene transmission spectra. 

The rock samples were observed at ambient temperature outside the 
laboratory in order to simulate as closely as possible practical remote sensing 
conditions. As these samples were mostly freshly broken rock, their reflectance 
can be much greater than the rocks observed in the airborne mode. Thus the 
contribution of a reflected sky spectrum in the data can not be ignored if 
similar measurements are to be made with the object of determining mineral content 
The data reduction procedure used in this analysis is a simplified form 
of that discussed for the airborne spectra. Thus the negligible path length 
enables one to make the approximation N sur » N sp . In the rest of the calculation 
the method is identical. 

4. Generation of System Responses Simulating MS PS Scanner Outputs 

The voltage (E) produced by any selected filter-detector combination in 
the infrared scanner is dependent upon the spectral response of this system (<KA)) 
and upon the energy incident upon the filter and is given by 

E /AX N sp (X) d X 04) 

In this analysis the spectral responses of the MSC multi-spectral scanner (MSDS) 
and the University of Michigan Scanner (see Figure A2) 
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were digitized at the same wavelengths at which the spectra over Pisgah Crater 
had been recorded, and the responses determined as if the MSDS were flying rather 
than the IR pallet in MX108. They were computed using Equation 5, 

E = Hi N sp (Ai) (5) 

using the radiance calculated as that incident at the spectrometer in the flights 
over Pisgah Crater. These spectral responses contain brightness temperature 
information, that is, they vary with their true temperature and their emittance. 
They differ thus from the spectral emittance normally used which have been 
corrected for temperature effects. 

The system responses and the various ratios which can be formed from 
them are given in the appendix. 



5. BMD07M Analysis of Data-MX108 


The BMD07M stepwise discriminant analysis program was used to. test 
several types of data derived from the airborne radiance measurements of MX108. 
The following transforms were used as "variables", 

1. 50 centermost wavelengths , in the spectral emittance 
curves . 

2. System Responses (9) , from the three channels used by 
Vincent and the 6 from the MSDS Scanner. 

3. Ratios (18), from 

(a) Vincent's channels (3), and 

(b) MSDS channels (15) 

(To reduce the calculations cross ratios between channels in (a) and (b) 
systems were not made. These, of course, can be theoretical ratios only.) 

Table A1 shows these variables in detail. 

Table A2 indicates the sequential success using the BMD07M analysis 
on the two "end-member" terrain types - mixed terrains (dry lake sediments, 
alluvium and a lava) and "all lava" terrains. The success, by rock type, is 
listed in detail in Appendix Table D5. 

Table A3 examines the sequence of variables chosen by the program, for 
analysis of each terrain, and by "method". The ranking in the initial or zero- 
step of BMD07M is quoted, that is, the "redundant" data is still included. When 
BMD07M selects its variable for the second and successive step this redundancy is 
minimized. It is unclear yet whether this is a desirable or undesirable concept. 

Figures A3 - A5 show the discriminant plots from the BMD07M analysis 
for the mixed terrain; A6 - A8 are for the all lava terrain. 



TABLE A1 


VARIABLES USED IN BMD07M ANALYSIS 


A. "50 Central Wavelengths" 


#1 = 
#50 = 

System Responses 

8.100 ym 
11.970 ym 



1 : 

Vincent #1 (VI) 

6 : 

19 


2 : 

V2 

7 : 

20 


3 : 

V3 

8 : 

21 


4 : 

5 : 

MSDS #17 

18 

9 : 

22 


Ratios 

1 : 

V1/V2 

10 : 


18/20 

2 : 

V1/V3 

11 : 


18/21 

3 : 

V2/V3 

12 : 


18/22 

4 : 

MSDS 17/18 

13 : 


19/20 

5 : 

17/19 

14 : 


19/21 

6 : 

17/20 

15 : 


19/22 

7 : 

17/21 

16 : 


20/21 

8 : 

17/22 

17 s 


20/22 

9 : 

18/19 

18. : 


21/22 
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TABLE A2 


PERCENTAGE CORRECTLY IDENTIFIED IN TRAINING SET (BMD07M) 


Mixed Terrain 
Dry Lake Seds. A (53) 

Alluvium A (27) 

P-Train Lava III (31) 

111 spectra 


All Lava Terrain 


Lava Phase I* ** 
Phase II 
Phase III 


(35) 

(62) 

(65) 

162 spectra 



A. Normal Method (Spectral Emittance) 


1. 50 Central A's 66 

B. Systems Responses 
1. Vincent's - 3 90 

2 ; , MSDS - 6 90 

r Ratios of Channel Outputs 

1. Vincent's - 3 86 

2, MSDS - 15 58 


*Phase I Lava 


(IE + IK) 


Phase II Lava ■ (IIA + IIB + IIC + IID + IIF) 


Phase III Lava 


(IIIA + III-TRAIN + IIIG + IIIH) 


**Success can oscillate with increased stepping 





TABLE A3 


INITIAL RANKING OF FIRST FOUR VARIABLES CHOSEN 
( Zero-Step of BMD07M ) 



Mixed Terrain 

All Lava Terrain 

A. EMITTANCE SPECTRA 

SEQUENCE 

SEQUENCE 

50 Central A's 

8.96; 8.80; 8.88; 9.04 pm 

11.75; 11.83; 11.07; 10.' 

B. SYSTEM RESPONSE 



Vincent's - 3 

VI; V3; V2; 

(VI = V2 = V3) 

MSDS - 6 

18; 17; 22; 21 

17; 20; (21 = 22) 


C. RATIOS 

Vincent's - 3 
MSDS - 15 


VI /V2; VI /V3; V2/V3; -- 
18/19; 18/20; 18/21; 17/21 


(VI /V2 = VI /V3) V2/V3 — 

17/18; 17/22; 17/21; 17/19 



23 0 -2.250 0.750 3.750 6.750 

-3.750 -0.750 2.250 5.250 8.250 
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BMD07M 

A. Summary - Effect of Temperature of the Target 

1. Temperature effects have been removed in spectral emittance calculations, 

wherein each spectrum is directly ratioed to a blackbody at its "apparent 
temperature" = 1.0 at one wavelength). 

2. Temperature effects are also minimized in the ratio calculations, 
wherein system responses of pairs of filter channels are ratioed. 

3. Original temperature effects still remain in the raw system response data. 

4. Because of these facts, the second most powerful discriminant is the 
rock radiance level ( system response data ) , which varies directly with the 
temperature of the target (alluvials and dry lake sediments are cooler than 
the black lavas). 

B. Summary - Spectral Emittance Data Only 

Using spectral emittance data (50 central wavelengths) temperature 
effects are minimal and chemical (raststrahlen) effects can dominate. Many of 
the closely-adjacent wavelengths are redundant (high correlation) coefficient in 
BMD, and hence not utilized, after the first step, by BMD07M. 

a. Mixed Terrain : A high success (> 90%) is achieved only after several steps. 
(See Figure A9 and Table D5) . Firstly Alluvium (A), and then Dry Lake Sediments (D) 
are correctly differentiated, but the P-Lavas take at least 4 steps to reach 75%. 

b. All Lavas : After the first steps Phase I (oldest) and Phase III (youngest) 
are classified at 50-60% level. Phase II reaches this level only after 4 steps. 

C. Summary - System Response (Radiances) Only 

a. Mixed Terrain : Responses from either Vincent's or MSDS channels reach 

90% correct- levels after 2 steps, much more quickly than the spectral emittance, 
hence are more simply diagnostic. 

b. All Lavas i Vincent's channel radiance levels are only 70-80% correct 
regardless of the number used. The MSDS system does better with Phase I and II 



iili 



Spectral Emittance Vincent's MSDS Vincent's MSDS 

Spectral Responses Ratios 
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lavas with increasing steps, but Phase III remains constant. 

D. Summary - Ratios 

a. Mixed Terrain : Vincent's ratios are correct at the 75-80% level at Step 1, 

but the MSDS system only reaches 77% after 2 steps. 

b. All Lavas : No clear pattern emerges; if anything extra steps decrease 

the success achieved. 

E. Conclusions 

Table A4 acts as a summary of these results. From this it is reasonably clear 
that the chemical effect sought is always effective as a discriminant in the first 
step of the classification process, but the success level varies with the method 
(system response - best; then spectral emittance, closely followed by the ratios). 
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TABLE A4 

INFORMATION USED IN DECISION SEQUENCE 


' ■ 

(BMD07M) 




* 

1 

MIXED TERRAIN 

ALL LAVAS 

Spectral Emittance 

Chemical 

? Temperature 

System Response 

1 


Vincent's 

Chem. -*■ Temp. -*■ Chem. 

Chem. - Chem. - Temp. 

MSDS 

Chem. -+■ Temp. -*■ Chem. 

Chem. - Temp. - Chem. 

Ratios 
Vincent's 
/ .. MSDS 

2k: 

. 

Chem. Temp . 

Temp. Chem. 

Chem. Temp. 

Chem. -*■ Temp. 














B. TASK 2.3 - 


Comparison of IR Spectral Emittance Data with K-Band Scatterometer Data . 

1. General Background 

The presentation of non- imaged data to an audience of users, oriented 
towards imagery, has always been a problem. Numeric data and data representing 
the changing values of one variable measured against time (or distance) form an 
increasing proportion of that collected by aircraft and satellites. As we 
progress towards a more clear understanding of the processes involved in the 
interaction of electromagnetic radiation with the earth, and towards predictive 
modelling, the need for numeric data rather than spatially displayed imagery 
becomes more evident. 

Attendant with the increase in numeric data is the problem of display, 
and the education of the audience to the manner in which significant changes are 
identified. One of the most difficult types of data to display are spectra. 

Yet these are capable of being gathered by an instrument like the infrared 
spectrometer at rates like six spectra per second. 

The purpose of this part of our program was to illustrate how we have 

adapted a colorizing method (developed for radar back-scatter spectra) to our 
infrared spectra, and show the comparison between radar spectra and infrared 
spectra obtained along the same flight line, albeit at different times. 

Both sets of data can be colorized in a comparable way as easy guides 

to visualizing their differing information contents. We have followed the lead 
of Moore et al. (1968^<*herein they displayed radar back-scattering spectra 
(a vs. look angle) as colors generated in a 3-gun color TV tube. 

(1) Moore, R.K. , Waite, W.P., Lundien, J.R. , and Masenthius, H.W. (1968) "Radar 
Scatterometer Data Analysis Techniques". Proc. 5th Symp. Remote Sens, of 
Environ., Ann Arbor, Mich., 1967, — (U of Michigan Press, Ann Arbor). 
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2. University of Kansas Presentation of Scatterometer Spectra 

Moore et al. (1968) found that these Scatterometer spectra (a vs 0) 

(Fig. Bl) could be rendered more meaningful to outside observers if a transform 
was performed involving subtraction of the curve from the group mean (Fig. B2) . 

A further and far more successful presentation was made after a Principal 
Components analysis revealed that 90% of the spectral variability would be 
accounted for by using only 3- of the 8- data points in the curve. These 3 
values were found to be the a value at 5 8 , 30° and 60° look angles (measured 
from the vertical, nadir = 0°). 

The values were converted numerically to equivalent voltages and used 
to produce colors in the color TV tube (Fig. B3) . 

3. Stanford Technique for IR Spectra 

a. Selection of 3-values using BMD07M 

The Stanford IR Spectra are composed of a large number of digitized 
emittance values. Our analysis using the Stepwise Discriminant program BMD07M 
enabled the 1-, 2- and 3- most significant values to be identified from this 
contribution to the discrimination process. These key wavelengths (X,Y and Z) 
are those at which the several rock types of the Pisgah area can be best 
separated. The program also identifies other wavelengths, which are "almost- 
as-significant", but only selects the most-significant pair or triad in the 
first 3 steps. A high degree of redundancy exists in the data and actually one 
or other of several closely-adjacent points could have been chosen with almost 

equal effect. 

2 

b. I S Digital Image Processes (MCFV) 

2 

The processer unit in the International Imaging System (IS) MCFV equipment 
has three sets of matrix boards which enable each color band to be created by 
16- levels each in the Red, Green and Blue color-guns in the nearby color TV tube. 







SIGMA VS THETA TWO 
FREQUENCIES 



COMBINATIONS SHOWN AS 
COLOR, LENGTH, WIDTH, AND 
ANGLE 


Figure B 3. Possible 
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MULTIPLE PASS CONTOURS 
RELATED BY COLOR 



ANGLE PROPERTY CONTOURS 
RELATED TO COLOR PROPERTY 
CONTOURS 


Modes of Presentation of Scatterometer Data. 
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Sixteen color bands can be created simultaneously, as parallel vertical bands 
in the tube. (This necessitates setting 16 trios of slides in the matrix board, 
each into one of 16 levels of intensity). (Fig. B4) . 

In this mode only the digital processing position of the unit is 
operated, and the front-end (TV camera and image processer) units turned off. 
Only the color display is activated and controlled by the digital settings on 
the matrix board panel. Some analog bias can be introduced to further saturate 
the colors, if desired. 

Each value of the spectral curve at the 3 (X,Y,Z) wavelengths is obtained 
and reconverted into a 16-level (4 bit) number. In our examples the spectral 
curves have been normalized (mean, set to 0.0 and S.D. = 1.0) so the values of 
X, Y and Z ranged from -2.0 to +2.0 before the transform and from 0 to 16 
afterwards. For each spectrum three new values of X*, Y* and Z* were calculated 
by this process. 

c. Operating Matrix Boards 

Almost infinite flexibility exists in the manner of assigning the X*, Y* 
and Z* values to the red, green and blue color guns. After considerable 
experimentation the following assignments were found most useful. 


CVF 

TABLE B-l 

ASSIGNMENTS 

Specific 

Counter Pt. 

Wavelength 

Descriptor 

Color-Gun Activated 

166 

8.565 ym 

X 

Blue } 

} 

Green JSlides 2, 3, 5 

121 

8.965 

Y 

153 

11.525 

Z 

Red } and 6 

or 128 

9.53 

Y' 

Red (special case used 
in Slide 9) 


Thus the silicic and felsic rocks with lower spectral emittance in 8-9.5 lim 
regions would show less blue and more red and green; while the femic basalts, with 
higher emittances at X and Y, would show more bluish, with less red. 



MCFV DIGITAL IMAGE PROCESSOR 


Fig. B4. Digital Processor: The very high speed (4GMHz) Digital 

Processor accepts the outputs of the A/D Converter. A transformation 
is set up in a matrix, giving each individual level slice a specific 
combination of red, green and blue signal weighting. Thus, for each 
channel, the 16 levels are each assigned a red, green and blue value. 

The color assignment is made by a switch with a discrete number of 
settings. These settings are infinitely repeatable. The output of the 
color matrix is decoded to a 4-bit code and converted to an analog signal, 
where it is scaled to give an independent gain control for that color in 
that channel. The switch is the point where the matrix conversion is 
actually performed. Each of the one-of-sixteen decoders feeds three 16 
by 16 matrix switches, one for each color. Since only one line is active 
at any time, the same output may be tied to a number of inputs, but the 
converse is not possible. Any level may be transformed in this manner to 
give any color, and, in fact, all levels in one channel may be transformed 
to the same color. The switch outputs are decoded and scaled as described 


above . 
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It should be noted that these assignments are completely arbitrary and 
any other observer could have created a differing (but parallel) set of false 
colors by choosing a different X-*, Y-*, and Z-*to-color gun associates. 

Upon setting the red-green-blue values for each spectrum a vertical band 
of color tone appeared on the TV tube representing the relative intensities of 
these 3 spectral emittance values. Color photographs (slides) were taken of 
each set of 14-16 bands so programmed, and used for data recording. (Usually 
only 14 were used to a black (non-activated) border on each side) . 

Several sets of bands have been color printed and are shown as Figures 
B5-9. Unfortunately, even today color prints do not have the full brilliance 
and saturation of the slides or the original screen, hut they give an idea of 
the color discrimination which can be produced by this method. 

d. Application to Pisgah Crater, California flight data (MX108) ♦ 

Fig. BIO shows the original color page from Moore et al. (1968) paper, 
on which have been placed the chips of the color prints prepared from the slides 
for set 5. This set of 14 spectra represents the most complete sampling of rock 
types possible during the flight line. 

e. Scatterometer Ground Pattern or "Footprint" 

The 2-cm, vertically-polarized scatterometer operating in the NASA 
aircraft on Mission 21 (MX21) has a fan-shaped beam 2.5° wide by 120° (fore and 
aft), pointing normal to the ground surface (Fig. Bll) . Discrete areas from 
which a spectrum of back-scatter (a) values may be obtained as a function of 
look-angle (from +60° to -60°) . The cross-track width is set by the 2° fan and 
is 140 feet at the 4000-4200 feet flight altitude used in Line 1, Run 3, Flight 5 
of MX21. The long- track dimension is a function of the post-flight Doppler 
filtering of the returned signals, and was set to give a square cell at the 
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+30° (forward view from the nadir) . This increases the square cell (or"foot- 
print")of the data block to 165 x 165 feet at the +30° viewing angle (Fig. B12) . 

The flight line across the lava flows and adjoining sedimentary and 
lake beds is shown in Fig. BIO (marked "Sc"). The line contained 135 of these 
footprints and was therefore further subdivided on a geological basis into the 
14 units shown therein. 

f. Infrared Spectrometer Ground Patterns 

In a comparable manner the 6-12 ym spectrometer operating from the 
NASA aircraft on Mission 108 (MX108) collected data from a 0.4° (7mrad) circular 
field of view, swept forward by the motion of the aircraft down the two flight 
lines shown also on Fig. B10 (marked IR 1701 and 1702), almost exactly parallel 
to that of the scatterometer. In a desert environment like this little change 
in terrain would be expected between the MX flight of 1966 and that of the MX108 
flight in 1968, except perhaps in the lake mud, which would change in texture 
each wet season. 

The field of view or ground pattern thus at the 2000 feet flight altitude 
would be 15 feet wide. Each spectrum takes 150msec to collect, during which the 
aircraft has moved forward about 40 feet, or 40 x 15 feet of terrain for each 
spectrum. The group mean spectra tabulated in Table B-2 are shown in Figure B13, 

Again the spectra were segregated on a geological basis into the groups 
shown by the 31 wider black bars along the two flight lines (Fig. BIO). The 
numbering has been modified so that the circled numbers are now directly 
relatable to those used in the scatterometer analysis. 

Generally 15-30 spectra, covering 15 feet width x 600-1200 feet length, 
were used in each group. In a similar pattern to the scatterometer path, the 
IR path often crossed changing lithologies before this optimal number of spectra 
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TABLE B2 


DESCRIPTIONS OF SPECTRAL GROUPS - ON FLIGHT 1 


MX108-1-PISGAH 
LOC* NAME 


NO. OF 

GMT 

GMT 

GROUP 

SPECTRA 

START 

STOP 

NUMBER 


la. 

Alluvium C 

25 

18:49 

44078 

18:49:51360 

1 

lb. 

Alluvium AC 

30 

18:50 

22915 

18:50:31410 

2 

2. 

Sand over Basalt II-C 

11 

19: 9 

23685 

19: 9:26719 

3 

3a. 

Alluvium A 

26 

18:49 

58035 

18:50:05923 

4 

3b. 

Alluvium B 

30 

18:50 

35050 

18:50:43864 

5 

4b. 

Sand over Basalt I-B 

21 

19: 8 

37569 

19: 8:43636 

6 

4a. 

Sand over Basalt I-A 

9 

18:50 

52041 

18:50:54468 

7 

6. 

Pisgah Flow III 

6 

19: 8 

27252 

19: 8:28784 

8 

7. 

Pisgah Flow II 

4 

19: 8 

25750 

19: 8:26647 

9 

8 . 

Pisgah Flow I 

15 

19: 8 

20896 

19: 8:25143 

10 

9a. 

Lava Flow II-A 

13 

18:51 

04785 

18:51:13887 

11 

9d. 

Lava Flow II-D 

9 

19: 8 

12083 

19: 8:16649 

12 

10a. 

Pisgah Lava III-A 

8 

18:51 

16921 

18:51:22383 

13 

lOtr. 

P-Train Lava (not included) 

31 

19: 7 

57521 

19: 8:06623 

14 

10b. 

Pisgah Lava III-B 

8 

18:51 

23596 

18:51:28452 

15 

10c. 

Pisgah Lava III-C 

18 

18:51 

29058 

18:51:34231 

16 

10g. 

Lava III-G 

10 

19: 7 

48418 

19: 7:55701 

17 

lOd. 

Pisgah Lava III-D 

12 

18:51 

34837 

18:51:38161 

18 

10h., 

Lava III-H 

8 

19: 7 

40241 

19: 7:46309 

19 

I6f . 

Pisgah Lava II-F 

20 

18:51 

47263 

18:51:53332 

20 

10 j. 

Lava III-J 

16 

19: 7 

31746 

19: 7:39317 

21 

I5d. 

Alluvium D 

21 

18:51 

58186 

18:52:04254 

22 

11. 

Lava I-K 

16 

19: 7 

15652 

19: 7:20506 

23 

16c. 

Lava II-C 

10 

19: 7 

: 08371 

19: 7:14439 

24 

16b. 

Lava Flow II-B 

10 

18:52 

12461 

18:52:16998 

25 

12. 

Pisgah Lava I-E 

19 

18:52 

19136 

18:52:25204 

26 

14b. 

Dry Lake Sediments B 

23 

19: 6 

44098 

19: 6:50774 

27 

14 a. 

Dry Lake Sediments A 

53 

18:57 

38314 

18:57:54090 

28 

14c. 

Dry Lake Sediments C 

25 

19: 6 

34390 

19: 6:41671 

29 

15f . 

Alluvium F 

42 

19: 6 

14078 

19: 6:26501 

30 

15e. 

Alluvium E 

27 

18:52 

59793 

18:53:07681 

31 


♦Locations on Fig. B 14 




Fig. B13. Nine sets of normalized 
group mean spectral emittance curves, 
all with the same horizontal scale 
(wavelength from 8.0 to 12.0ym). The 
vertical scale in all sets is constant, 
but the sets have been displaced vert- 
ically for ease in viewing. Normalized 
spectral emittance is plotted as the 
ordinate, with the mean ■ 0.0 as the 
central line on each spectrum. 

Spectral shapes thus may be directly 
compared. 




accumulated (especially near the crater) . The group then contains less spectra 
(N, varying from 4-15) with correspondingly higher relative error in the 
estimation of the group mean spectrum. 

TABLE B-3 DIMENSIONS OF SPECTRA ON GROUND 



MX21 (4/5/66) 

MX108 (10/8/68) 

1 

Scatterometer 

1 

IR Spectrometer 


F.O.V : 

2° wide 

0.4° circular 


Flight Alt. Terrain : 

4000-4200 ft. one flight 

2000-2500 ft. two lines 



line (103) 

(1701, 1702) 


Ground patch width : 

165' (+30° view) 

15' (vertical) 


Ground patch length : 

165' 

40' 


// of Spectra : 

i 

168 Spectra/14 groups 

600 Spectra/31 groups 


4. Results 


a. The method of colorizing appears to emphasize differences in the 
spectral curves quite clearly. 

b. It is simple to achieve providing one has access to the specific units 
(I 2 S - MCVF or the University of Kansas equivalent) . 

c. It is not clear however if the method can be developed (without a 
further research period) to give suitably diagnostic colors of a unique tone for 
each setting without additional 3 color— saturation from the analog processer. 
Above each color slide print we have indicated the analog settings for the red— 
blue-green guns as for example in Slide 9 (R[10x] ; G[4.7x]; Bl2.8xJ). It will be 
noted that those vary from slide to slide and throw a non— linear bias into our 
comparisons. 

d. The 14 color pattern shown in Rig. BIO and tabulated (Tahle B-4) do 
show the same relationship between the scatterometer color sequence and the IR 
color sequence, indicating a close general correlation between the two methods. 









TABLE B- 5 SUMMARY 
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COLORS PRODUCED 


SCATTEROMETER 

INFRARED 

SIMILARITIES 



1. Sands and Alluvials 

Browns and Tan 

Red, Orange and 

(#1-3, 14, 15) 


Brown 

2. Basalt Flows 

Blue, Grey-Blue and 

Blue, Light Green 

(#5-12) 

Grey 

Light Blues 

DISSIMILARITIES 



1. Sands over Basalt 

Grey Tan 

Grey-Blue-Green 

(#4a, 4b) 


Avocado Green 

(#2) 

Dark Tan 

Red 

2. Lake Beds 

Dark Tan 

Very Light Blue 

(#14) 

3. Specific Lava Flows 



(#7) 

Grey 

Fuchsia 


e. The dissimilarities are indicative of the effect of surface roughness 
(scatterometer) and the shallow "skin depth" or penetration shown by the IR data. 
The windblown sands (Site 2. and Sites 4a and 4b) are cases in point, showing 
more of the character of sands to both the IR and scatterometer. With the 
various lava flows and cinders of the Pisgah Crater area the scatterometer 
generally agrees with the IR (blues and greys) . Some lava flows C#7) are 
different in the IR data, being more like alluvlals. 

The most clear dissimilarity is in the dry lake sediments which the 
scatterometer sees more as alluvlals (light grey tan and dark tan) while the IR 
data closely agree with the basalt composition. (This has been further 
substantiated by the work of Vincent and Thompson using emittance-ratio imagery) . 



5. Conclusion 

The method is a valid and dramatic technique for displaying spectral 
data. It needs further study if one wishes to make unique colors from 
digital input alone. 



57 


C. Standard Infrared Spectral File 
1. MX108 Airborne Spectra 

The following spectra were analyzed as described in Section A3 (a) . The 
spectra are plotted between 8 and 12 microms as it was felt that this is the 
only region where atmospheric effects could be reliably eliminated. The scale 
is such that each step is 0.001 of an emittance unit. The heavy curve is the 
average "effective emissivity" for the group in question. The other curve is 
the standard deviation of this main curve. The table at the top gives the number 
of spectra in the group, the location and time of measurement and the average 
temperature with its standard deviation as calculated from the spectral data. 

It should be noted that full scale on the graph is not equal to 1 but 
that the effective emittance and the standard deviations are plotted from 1 and 
0 respectively. 


f 



TABLE Cl 



DESCRIPTIONS OF 

SPECTRAL GROUPS - ON FLIGHT 

1 LINE 




AIRBORNE 

SPECTRA 




MX108-1-PISGAH 







NO. OF 

GMT 

GMT 

LOC 

NAME 

SPECTRA 

START 

STOP 

1 . 

Alluvium C 

25 

18:49:44078 

18:49 

51360 

2. 

Alluvium AC 

30 

18:50:22915 

18:50 

31410 

3. 

Sand over Basalt II-C 

11 

19: 9:23685 

19: 9 

26719 

4. 

Alluvium A 

26 

18:49:58035 

18:50 

05923 

5. 

Alluvium B 

30 

18:50:35050 

18:50 

43864 

6. 

Sand over Basalt I-B 

21 

19: 8:37569 

19: 8 

43636 

7. 

Sand over Basalt I-A 

9 

18:50:52041 

18:50 

54468 

8. 

Pisgah Flow III 

6 

19: 8:27252 

19: 8 

28784 

9. 

Pisgah Flow II 

4 

19: 8:25750 

19: 8 

26647 

10. 

Pisgah Flow I 

15 

19: 8:20896 

19: 8 

25143 

11. 

Lava Flow II-A 

13 

18:51:04785 

18:51 

13887 

12. 

Lava Flow II-D 

9 

19: 8:12083 

19: 8 

16649 

13. 

Pisgah Lava III-A 

8 

18:51:16921 

18:51 

22383 

14. 

P- Train Lava (not included) 

31 

19: 7:57521 

19: 8 

06623 

15. 

Pisgah Lava III-B 

8 

18:51:23596 

18:51 

28452 

16. 

Pisgah Lava III-C 

18 

18:51:29058 

18:51 

34231 

17. 

Lava III-G 

10 

19: 7:48418 

19: 7 

55701 

18. 

Pisgah Lava III-D 

12 

18:51:34837 

18:51 

38161 

19. 

Lava III-H 

8 

19: 7:40241 

19: 7 

46309 

20. 

Pisgah Lava II-F 

20 

18:51:47263 

18:51 

53332 

21. 

Lava III-J 

16 

19: 7:31746 

19: 7 

39317 

22. 

Alluvium D 

21 

18:51:58186 

18:52 

04254 

23. 

Lava I-K 

16 

19: 7:15652 

19: 7 

20506 

24. 

Lava II-C 

10 

19: 7:08371 

19: 7 

14439 

25. 

Lava Flow II-B 

10 

18:52 12461 

18:52 

16998 

26. 

Pisgah Lava I-E 

19 

18:52:19136 

18:52 

25204 

27. 

Dry Lake Sediments B 

23 

19: 6:44098 

19: 6 

50774 

28. 

Dry Lake Sediments A 

53 

18:57:38314 

18:57 

54090 

29. 

Dry Lake Sediments C 

25 

19: 6:34390 

19: 6 

41671 

30. 

Alluvium F 

42 

19: 6:14078 

19: 6 

26501 

31. 

Alluvium E 

27 

18:52:59793 

18:53 

07681 

44. 

Pisgah Cinders I 

4 

18:50:55076 

18:50: 

56000 

45. 

Pisgah Cinders II 

5 

18:50:56607 

18:50: 

57821 

46. 

Pisgah Cinders III 

9 

18:50:58427 

18:50:00855 

54. 

Pisgah Lava III (A-D) 

53 

18:51:21169 

18:51: 

36948 


MX108- 1-SUNSHINE 





40. 

Sunshine Lava A 

22 

18:57:11905 

18:57: 

18290 

41. 

Sunshine Lava B 

25 

18:57:04017 

18:57: 

11297 

42. 

Sunshine Cinders C 

17 

18:56:54626 

18:56: 

59480 

43. 

Sunshine Cinders D 

20 

18:56:47633 

18:56:53412 


VARIOUS 





53. 

Palmdale Lake 

44 

17:12:38496 

17:12: 

52143 


53 
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C2. Ground Spectra - Short Path Length 

These spectra were analyzed as described in Section A3(b). The basic 
format is the same as for the airborne spectra; however the spectra have been 
digitized at regular wavelength intervals (0.05y) providing more points. The 
standard deviations have not been included as the extremely low noise level 
provided by the ground based equipment in comparison with airborne spectra over 

a fluctuating terrain would not give a meaningful comparison. 

The following table contains a brief mineralogical description of the 

rocks studied. Not all the rocks described are included in the spectra due to 
experimental difficulties at the time of measurement. The spectra can be 
correlated with the table through the time indicated at the top of each spectrum. 
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TABLE C2 

DESCRIPTION OF SAMPLES USED FOR GROUND SPECTRA (Site #27) 


Run Sample Time 


Mineralogy of the 1 1/2" x 1 1/2 11 Sample Area Surface 


1. 10:40 

Cinko Lake 
Granodiorite 


The surface is coated approximately 75% by fine Rough 
grained, black tourmaline crystals. The remain- 
ing lighter area is largely quartz 20% and feld- 
spar 5%. 


2. 10:45 

Cinko Lake 
Granodiorite 


3. 10:50 

Cinko Lake 
Granodiorite 


The sample area is approximately 30% quartz — Polished 
30% biotite, 20% feldspar mostly plagioclase, 

15% hornblende 5% accessory minerals. The 
texture is medium grained equigranular-grano- 
diorite. 

The 1/2 M xenolith in a matrix of Cinko Lake Sawed 

granodiorite is composed of fine grained biotite 
80% and hornblende 20%. 


4. 11:00 The sample is approximately the same composition Sawed 

Cinko Lake as that seen in Run #2 at 10:45. 

Granodiorite 


5. 11:05 

Cinko Lake 
Granodiorite 

6 . 11:10 


Fremont Lake 
Granodiorite 


The surface is coated approximately 50% by fine Rough 
grained, dark tourmaline crystals. The lighter 
material is approximately 25% quartz and 25% 
feldspar. 

The sample area is medium grained, equi granular, Rough 
and is composed of 60% feldspar, mostly plagio- 
clase, 15% quartz, 10% biotite, 10% hornblende, 
and 5% accessory minerals. 


7. NASA #302 11:15 
Cascade Creek 
Granite 


8. NASA #162 11:25 
Dorothy Lake 
Alaskite-Granite 


The sample area is medium grained hypodiomorphic. Rough 
and is composed of 50% quartz, 30% orthoclase, 

10% biotite, 2% hornblende, some plagioclase and 
accessory minerals. 

The sample is a fine grained texture composed of Sawed 
60% feldspar, 30% quartz, 2% biotite and accessory 
minerals . 


9. NASA #162 11:30 The sample is approximately the same as Run #8 Rough 

Dorothy Lake at 11:25. 

Alaskite-Granite 


10. NASA #308 11:35 
Millcreek 

Porphyritic Quartz Mon- 
zonite 


The sample has porphyritic phenocryst of ortho- Rough 
clase in a coarse grained matrix of 30% ortho- 
clase, 20% plagioclase, 30% quartz, 10% biotite, 

5% hornblende and 5% accessory minerals. 
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Run Sample Tine 

11. NASA if 308 11:40 
Millcreek 

Porphyritic Quartz Mon- 
zonite 


Mineralogy of the 1 1/2" x 1 1/2" Sample Area Surface 


The xenolith is composed of approximately 50% 
biotite and 50% quartz in a very fine grained 
equigranular matrix. 


Rough 


12. NASA if 316 11:45 
Patterson Grade 
Granodiorite 


The sample area is medium grained, equigranular Rough 
and composed of 60% feldspar mostly plagioclase, 

15% quartz, 15% hornblende and biotite and 10% 
accessory minerals. 


13. NASA if 383 11:50 
Cathedral Peak 
Porphyritic Quartz Mon- 
zonite 


The sample has orthoclase (?) phenocryst in a 
coarse grained matrix of feldspar approximately 
30% orthoclase and 30% plagioclase, 30% quartz, 
5% biotite and 5% accessory minerals. The 
surface area was moderately weathered. 


Rough 


14. NASA if 383 11:55 


Same as above, except fresh rathered than Rough 

weathered. 


15. NASA #331 12:05 
Topaz Lake 


1 1/2" microcline phenocryst in a matrix of Run Rough 
#17, Sample #331 at 12:20. 


16. NASA if 331 12:15 
Topaz Lake 


1 1/2" microcline phenocryst in a matrix of the Rough 
below sample - NASA #331. 


17. NASA if 331 12:20 
Topaz Lake 

Porphyritic Quartz Mon- 
zonite 

(General pass) 


The sample has microcline phenocryst in a matrix Rough 
of coarse grained subhedral crystals composed of 
35% microcline, 30% plagioclase, 25% quartz, 4% 
biotite, 6% accessory minerals. 


18-19-20 Calibration 


21. NASA #621 14:50 
Brown Bear Pass 
Basalt 


Weathered surface of basalt, some hematite Rough 

staining. 


22. NASA #621 15:00 
Brown Bear Pass 
Basalt 


Fresh surface is composed of 70% plagioclase Rough 

feldspar, 15% augite, 5% orthoclase, 5% pyroxene 
and 5% magnetite weathering to hematite. 


23. Q #8 15:05 

Crow Springs 
Porphyritic Quartz Mon- 
zonite 


The sample has porphyritic phenocrysts of plagio- Rough 
clase 30%, and interstitial quartz 25%, ortho- 
clase 30% and hornblende 10% and 5% accessory 
minerals . 


24. Q #18 15:15 Sample has medium grained matrix, phenocrysts Sawed 

Crow Springs mostly well-formed plagioclase 35% up to 5mm. in 

Quartz Monzonite Porphyry in length, quartz 25%, orthoclase 30%, biotite 50%, 

some hornblende, and the rest accessory minerals 
(Dark Phase). 
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Run 

Sample 

Time 

25. 

Q #18 

15:25 

26. 

Q #50 

15:30 


Crow Springs 


27. 

Q #71 

15:35 


Crow Springs 

28. 

Q #71 

15:40 

29. 

Q #1 

15:45 


Crow Springs 

30. 

Q #61 

15:50 


Crow Springs 

31. 

Q #77 

16:05 


Crow Springs 

32. 

Q #77 

16:10 

33. 

Q #63 

16:15 


Crow Springs 

34. 

Q #63 

16:20 


Crow Springs 

35. 

Q #56 

16:25 


Crow Springs 

36. 

Q #58 

16:35 


Crow Springs 

37. 

Q #58 

16:40 

38. 

Q #70 

16:45 


Crow Springs 


Mineralogy of the 1 1/2" x 1 1/2" Sample Area 

Same as above Q #18 - 15:15 except it has a 
rough surface. 

Sample area strongly welded, ash flow tuff, 
completely devitrified, axiolitic texture, 
composition 60% glass and ash devitrified to 
cristobalite and K-feldspar. 15% subhedral 
quartz, 10% sanidine with trace of biotite and 
magnetite. 

Fine grained (hypocrystalline) with microlite 
matrix. Approximately 50% of area mostly plagio- 
clase, larger plagioclase, subhedral to euhedral 
(21%), augite 8%, glass 17%. 

Same as above except deeply weathered, magnetite 
is forming irons tain. 

Strongly welded quartz latite. Composition - 
30% plagioclase, 10% quartz, 10% biotite, in a 
matrix of 50% devitrified glass. 

Non-welded lithic tuff. Composition - 50% 
volcanic dust, 14% subhedral sanidine, 12% 
quartz, 2% biotite, 10% pumice fragments. 

Weathered vitrophere, strongly welded, squashed 
fiamine filled with glass fragments. 

Non-weathered side of the above sample. 

Welded quartz latite - 20% plagioclase, 15% 
sanidine, 10% quartz, 5% biotite, some hornblende 
- 10% fiamine. The matrix is composed of 40% 
devitrified shards. The sample is weathered. 

Approximately the same as above except the 
sample is fresh rather than weathered. 

Strongly altered obsidian or welded tuff - 
strongly devitrified 40% glass, 25% cristobalite, 
20% sericite (?), 10% feldspar, 5% quartz. 

Strongly welded crystal tuff. 15% sanidine, 10% 
quartz, 10% fiamine. The matrix is composed of 
60% glass shards which have been devitrified. 

Same as above except for sawed surface. 

Strongly welded ash flow tuff 5% plagioclase, 

5% sanidine, 2% quartz, 20% lithic fragments, 

68% severely welded glass shards - reddish brown, 
devitrified to cristobalite and K-feldspar - a 
weathered sample. 
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Surface 

Rough 

Rough 

Rough 

Rough 

Rough 

Rough 

Rough 

Rough 

Rough 

Rough 

Rough 

Rough 
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Rough 
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0.440 o. j> i 0.>a/ 

• o.469 b. /'.a 

ti.'liiO 0.4*0 to. if. « J./**.* 

4*/t/U J.4*d V*2«o 

to » 4 ‘ * / 6 . I* to w j . # 0 to 

iu.jgj 0.934 UiJ'/.i J.'/tiU 

t ««•**«.,} J • 4 7 1 . J . */ / ,* 

I .«)!) 0.4/6 | U. >t j J J./// 

i I */wJ u.//d { 1 i / *J j o • 4 / 4 

II . w«ij u ll.o to J j • Jo j 

i^*o iu 0**» t»u 


to • 1 to 

L> * '/to'* 

ti. 1 63 

0.41*2 
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4oj 
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b - 4 6 j 

o .4 to 1 

/ ( i O to 

w. 'JJ O 

4 . * to J 
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4 a t\* w 

U.-/I.7 

4 . f 6 J 

o • 4 5 / 

t J. iJj 

u. 4to J 

i J. 1 too 

U.4t.2 

I U . to Jto 

0.//^ 

1 0. Ou 

0.4/4 

i J. i.UJ 

J . 4 7 / 

l 0 * 9 5 to 

0 . 9 / 7 

i 4 • jOto 

0.4/ / 

1 i . J 5 to 

0.9 /*> 

1 1. /oO 

J . 4 to 4 

11. /too 

0. 9 c* 6 


fj.2jj o.ojj u.250 0 . 4/3 

u. COO J.Vjto b.toOJ o.'/tou 

O.OjO 0.951 S.otoJ 0.961 

‘>•*00 O.Oto** O.'tOJ g.*5t> 

V. bOO o.4to4 f/.otoo 0.464 

il./vlJ 0 . 9l* 4 10.250 u . v O / 

1 J. Ljj (J. '//’* lO.totoJ 0.4/4 

1I.UC0 0.0// li.JbO 0.4/7 

11.400 J.to/V i 1 . 4 to J 0.474 
U.3JJ J. 9<*4 ll.oOO 0.4b* 


a. 000 0.972 a. 350 0.970 

b. /OO 0.955 8*750 0.956 
4. 100 0.449 9.150 0.9*6 
4.500 0.95b 9.55C 0.957 
4.403 0.959 9.950 0.958 

10.300 0.*6d 10.350 0.966 
*0./u0 0.974 10.750 0.975 
11. 100 0.977 11.150 0.978 
11.600 0.960 ll.550C.980 
11.900 0.985 11.950 0.985 
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t 3*1 
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I 8.3 
I 8.4 
1 8*4 
I 8.5 
I 8.5 
I 8.6 
f 8.6 
I 8. 7 
i 8.7 
I 8. 8 
I 6.8 
f 8.9 
i 8.9 
{ 9.0 
I 9 # 0 
I 9.1 
I 9. 1 
I 9.2 
I 9.2 
I 9.3 
I 9.3 
I 9.4 
I 9.4 
I 9.5 
I 9.5 
I 9.6 
I 9.6 
I 9.7 
I 9.7 
I 9.8 
I 9.8 
I 9.9 
I 9.9 
1 10. Q 
110.0 
lio.l 

1 10. I 
I 10.2 
l IQ .2 
I10.3< 
110.3 

1 10. 6( 
110.4: 
U0.5< 
110.5! 
|10.6< 
110.6! 
(10.7C 
110.7! 

I io. at 

110. 8f 
110.90 
I 10.95 

111. QC 
I 11.06 
I 11.10 
ai.15 
m.2o 

I 11.25 
111.30 
111.35 
111.40 
111.65 
111.60 
111.55 
111.60 
I 11.65 
I 11.70 
(11.75 
(11.80 
(11.85 
111.90 
(11.95 
(12.00 
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u . * 5 O 

0.9JJ 

6 • JOO 

0.920 
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* 
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(< • t< J J 
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0.411 

6. OL J 

b. 4 Jo 
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0.898 
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I 3.15 
I 8.20 
I 3.25 
I 8.30 
I 8.35 
i 8.40 
I 8*45 
I 8.50 
I 8.55 
| 8.60 
I 8.65 
I 8.70 
I 8.75 
1 8.80 
I 8.85 
1 8.90 
I 8.95 
1 9.00 
I 9.05 
| 9.10 
I 9.15 
I 9.20 
1 9.25 
I 9.30 
I 9.35 
I 9.40 
I 9.45 
J 9.50 
I 9.55 
1 9.60 
I 9.65 
I 9.70 
I 9.75 
I 9.80 
1 9.35 
I 9.90 
1 9.95 
110.00 
I 10.05 
110.10 
110.15 
110.20 
I 10.25 
110.30 
1 10.35 
110.40 
110.45 
I 10.50 
110.55 
110.60 
110.65 
I 10.70 
110.75 
I 10.80 
110.85 
110.90 
110.95 
111.00 
*11.05 
111*10 
1 11.15 

Hi. 20 

111.25 
111.30 
111.35 
(11.40 
1 11.45 
111.50 
I 11.55 
111.60 
111.65 
I ll* 70 
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I 11.80 
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111.90 
I 11.45 
| 12.00 
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APPENDIX TABLE Pi 
FILTER TRANSMISSION FUNCTIONS 

Wavelength Channels 



*17 

*18 

#19 

#20 

#21 

VI 

V2 

8.02 

— 

— 

— 

— 

— 

.06 

.01 

8.10 

.02 

— 

— 

— 

— 

.11 

.01 

8.18 

.18 

— 

— 

— 

— 

.18 

.01 

8.26 

.41 

— 

— 

— 

— 

.35 

.02 

8.32 

.57 

— 

— 

— 

— 

.43 

.02 

8.41 

.78 

— 

— 

— 

— 

.49 

.02 

8.48 

.90 

— 

— 

— 

— 

.53 

.02 

8.56 

.99 

.03 

— 

— 

— 

.62 

.02 

8.64 

.90 

.03 

— 

— 

— 

.69 

.02 

8.72 

.65 

.12 

— 

— 

— 

.72 

.02 

8.80 

.37 

.36 

— 

— 

— 

.76 

.02 

8.88 

.15 

.57 

— 

— 

— 

.77 

.02 

8.96 

.05 

.75 

— 

— 

— 

.78 

.02 

9.04 

.01 

.95 

.01 

— 

— 

.78 

.02 

9.13 

— 

.99 

.02 

— 

— ' 

.78 

.03 

9.20 

— 

.91 

.05 

— 

— 

.77 

.03 

9.27 

— 

.73 

.17 

.02 

— 

.77 

.04 

9.37 

— 

.42 

.42 

.03 

— 

.76 

.07 

9.45 

— 

.14 

.68 

.03 

— 

.74 

.12 

9.53 

— 

.03 

.96 

.04 

— 

.72 

.15 

9.65 

— 

— 

.92 

.04 

— 

.68 

.20 

9.69 

— 

— 

.83 

.05 

— 

.61 

.26 

9.77 

— 

— 

.64 

.06 

— 

.56 

.32 

9.86 

__ 



.42 

.11 

— 

.49 

.40 


D 1 


V3 



APPENDIX TABLE D1 (cont'd) 
FILTER TRANSMISSION FUNCTIONS 

Wavelength Channels 


10.01 

J 17 

#18 

#19 

.10 

1 20 
.41 

#21 

10.10 

— 

— 

— 

.62 

— 

10.17 

— 

— 

— 

.83 

— 

10.25 

— 

— 

— 

.91 

— 

10.34 

— 

— 

— 

.97 

— 

10.42 

— 

— 

— 

.98 

— 

10.49 

— 

— 

— 

.98 

— 

10.58 

— 

— 

— 

.99 

— 

10.66 

— 

— 

— 

.92 

— 

10.74 

— 

— 

— 

.83 

— 

10.81 

— 

— 

— 

.74 

CM 

o 

• 

10.90 

— 

— 

— 

.66 

.08 

10.97 

— 

— 

. — 

.45 

.29 

11.07 

— 

— 

— 

.18 

.59 

11.15 

— 

— 

— 

.05 

.76 

11.22 

— 

— 

— 

— 

.91 

11.30 

— 

— 

— 

— 

.98 

11.38 

— 

— 

— 

% 

.97 

11.44 

— 

— 

— 

— 

.92 

11.52 

— 

— 

— 

— 

.73 

11.60 

— 

— 

— 

— 

.51 

11.68 

— 

— 

— 

— 

.59 

11.75 

— 

— 

— 

— 

.53 

11.83 

— 

— 

— 

— 

.43 

11.89 

— 

— 

— 

— 

.38 


137. 
D lCcont'd) 


VI 

V 2 

V3 

.40 

.52 

— 

.37 

.58 

— 

.33 

.64 

— 

.30 

.69 

— 

.26 

.75 

— 

.23 

.81 

— 

o 

CM 

• 

.82 

— 

.17 

.82 

— 

.16 

.79 

— ' . 

.13 

.74 

— 

.11 

.68 

— 

.09 

.61 

— 

.07 

.54 

.01 

.06 

.48 

.20 

.05 

.42 

.30 

.04 

.38 

.60 

.04 

.34 

.14 

.04 

.31 

.24 

.04 

.27 

.38 

.03 

.24 

.60 

.03 

.21 

.74 

.03 

.19 

.68 

.03 

.17 

.65 

.02 

.16 

.64 

.02 

.13 

.65 



Wavelength 

#17 

#18 

APPENDIX TABLE D 1 
FILTER TRANSMISSION 

Channels 
#19 #20 

(cont'd) 

FUNCTIONS 

#21 

VI 

D 1 
V2 

138. 

(cont’d) 

V3 










11.97 

— 

— 

— 

— 

.26 

.02 

.12 

.68 

12.04 

— 

— 

— 

— 

.15 

.02 

.11 

.70 

12.12 

— 

— 

— 

— 

.07 

— 

— 

.75 

12.19 

— 

— 

— 

— 

.15 

— 

— 

.77 

12.26 

— 

— 

— 

— 

— 

r- 

— 

.74 

12.33 

— 

— 

— 

— 

— 

— 

— 

.68 

12.40 

— 

— 

— 

— 

— 

— 

— 

.65 

12.47 


— 

— 

— 

— 

— 

— . 

.61 

12.54 

— 

— 

— 

— 

— 

— 

— 

.57 

12.60 

— 

— 

— 

— 

— 

— 

— 

.55 

12.67 

— 

— 

— 



— 

— 

— 

.50 

12.74 

— 

— 

— 

— 

— 

1 

— 

.48 

12.81 

— 

— 

— 

— 

— 

— 

— 

.44 

12.89 

— 

’ 

— 

— 

— 

— 

— 

.42 

12.95 

— 

— 

— 

— 

— i 

— 

— 

.40 

13.02 

— . 

— 

— 

— 

— 

— ' 

— 

.39 

13.08 

— 

— 

— 

— 

— 

— 

— 

.37 

13.15 

— 

— 

— 

— 

— 

— 

— 

.36 

13.22 

— 

— 

— 

— 

— 

— 

— 

.35 

13.28 

— 

— 

— 

— 

— 

— 

— 

.34 




This table 

contains 

and 

APPENDIX TABLE D2 

the calculated System response!' for the 
University of Michigan Scanners 

MSDS 

139. 

D2. 

GROUP 1 

17 

18 

Channel 
19 20 

21 

22 

VI 

■ V2 

V3 

— 

— 

■ ' 



— 



AVERAGES 

19.58 

16.28 

14.84 

6.43 

6.60 

6.05 

12.35 

9.76 

11.59 

STD..DEVS. 

0.09 

0.07 

0.06 

0.03 

0.04 

0.04 

0.06 

0.04 

C.06 

S ERROR 

0.44 

0.41 

0.39 

0.47 

0.64 

0.74 

0.50 

0.39 

0.55 

GROUP 2 

AVERAGES 

20.09 

16.56 

15.04 

6. 64 

6.80 

6.19 

12.57 

9.90 

11.75 

STD. DEVS. 

0.11 

0.07 

0.06 

0.05 

0.05 

0.04 

0.06 

0.05 

0.06 

X ERROR 

0.55 

0.45 

0.42 

0.6 8 

0.66 

0.73 

0.50 

0.46 

0.52 

GROUP 3 
AVERAGES 

20.00 

16.73 

15.28 

6.57 

6.71 

6.18 

12.70 

10.05 

11.94 

STD. DEVS. 

0.11 

0.11 

0.09 

0.06 

0.04 

0.06 

0.09 

0.07 

0.07 

X ERROR 

0.57 

0.65 

0.59 

0. 86 

0.56 

0.91 

0.72 

0.73 

0.60 

GROUP A 
AVERAGES 

19.38 

16.14 

14.73 

6.36 

6.52 

6.00 

12.25 

9.69 

11.50 

STD. DEVS. 

0.08 

0.06 

0.07 

0.04 

0.04 

0.03 

0.06 

0.04 

0.07 

X ERROR 

0.41 

0.39 

0.46 

0.62 

0.57 

0.57 

0.48 

0.46 

0.57 

GROUP 5 
AVERAGES 

20.49 

16.76 

15.22 

6.81 

6.95 

6.29 

12.71 

10.01 

11.89 

STD. DEVS. 

0.22 

0.17 

0.13 

0.09 

0.08 

0.07 

0.14 

0.09 

0.10 

X ERROR 

1.09 

1.03 

0.83 

1.33 

1.09 

1.07 

1.12 

0.93 

0.87 

GROUP 6 
AVERAGES 

21.63 

17.66 

15.84 

7.18 

7.35 

6.63 

13.41 

10.47 

12.36 

STD. DEVS. 

0.48 

0.38 

0.37 

0.21 

0.18 

0.15 

0.29 

0.2 7 

0.29 

X ERROR 

2.24 

2.15 

2.34 

2.87 

2.49 

2.19 

2.13 

2.58 

2.33 

GROUP 7 
AVERAGES 

20.62 

16.81 

15.21 

6.84 

7.01 

6.33 

12.76 

10.01 

11.88 

STD. DEVS. 

0.34 

0.38 

0.39 

0.10 

0.11 

0.11 

0.30 

0.25 

0.31 

% ERROR 

1.66 

2.27 

2.54 

1.50 

1.51 

1.73 

2.34 

2.47 

2.64 

GROUP 8 
AVERAGES 

22.22 

17.90 

15.96 

7.45 

7.57 

6.80 

13.59 

10.56 

12.46 

STD. DEVS. 

1.63 

1.2C 

0.77 

0.54 

0.57 

0.52 

0.94 

0.60 

0.57 

X ERROR 

7.31 

6.73 

4.82 

7.25 

7.54 

7.60 

6.91 

5.69 

4.54 

GROUP 9 

AVERAGES 

23.20 

18.81 

16.60 

7.74 

7.89 

7.11 

14.26 

11.10 

12.90 

STD. DEVS. 

0.80 

0.48 

0.34 

0.27 

0.32 

0.25 

0.35 

0.26 

0.25 

* ERROR 

3.44 

2.56 

2.05 

3.43 

4.08 

3.52 

2.46 

2.35 

1.97 

GROUP LO 
AVERAGES 

21.21 

17.10 

15.22 

7. 11 

7.24 

6.46 

12.99 

10.08 

11.87 

STD. DEVS. 

0.79 

0.57 

0.46 

0.27 

0.28 

0.26 

0.44 

0.29 

0.38 

X ERROR 

3.71 

3.33 

3.05 

3.86 

3.89 

4.09 

3.42 

2.93 

3.23 

GROUP 11 
AVERAGES 

20.40 

16.52 

14.92 

6.82 

6.97 

6.23 

12.53 

9.83 

11.64 

STD. DEVS. 

0.54 

0.39 

0.20 

0. 18 

0.18 

0.21 

0.30 

0.18 

0.15 

X ERROR 

2.65 

2.35 

1.35 

2.62 

2.60 

3.30 

2.37 

1.85 

1.33 

GROUP 12 

A V ERAGES 

21 .36 

17.25 

15.56 

7.16 

7.32 

6.49 

13.08 

10.26 

12.14 

STD. DEVS. 

0.35 

0.22 

0.29 

0.15 

0.13 

0.10 

0.16 

0.16 

0.24 

X ERROR 

1.63 

1.25 

1.87 

2. 11 

1.81 

1.47 

1.23 

1.53 

1.99 


* Watts cm -2 ster -1 X10 3 
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140. 

D2 <Qont'd) 


GROUP 13 

17 

18 

19 

20 

21 

22 

VI 

V2 

V3 

AVERAGES 

20*29 

16.40 

14.79 

6.78 

6.92 

6.20 

12.45 

9.72 

11.57 

STO.DEVS, 

0.42 

0*30 

0.26 

0.17 

0.16 

0.11 

0.22 

0.18 

0.20 

* ERROR 

2.07 

1.80 

1.73 

2.44 

2.30 

1.81 

1.78 

1.87 

1.77 

GROUP 14 
AVERAGES 

21.29 

17*19 

15.52 

7.15 

7.26 

6.48 

13.05 

10.21 

12.13 

STO.DEVS, 

0.62 

0*43 

0.29 

0.26 

0.22 

0.21 

0.36 

0.22 

0.23 

X ERROR 

2.92 

2.53 

1.84 

3.68 

3.06 

3.19 

2.75 

2.15 

1.86 

GROUP 15 
AVERAGES 

20*48 

16.59 

14.96 

6.86 

6.98 

6.24 

12.60 

9.84 

11.69 

STO.DEVS. 

0.23 

0*18 

0.20 

0.10 

0.10 

0.06 

0.14 

0.12 

0.18 

X ERROR 

1.13 

1.11 

1.31 

1.41 

1.37 

1.04 

1.14 

1.21 

1.51 

GROUP 16 
AVERAGES 

20*22 

16.32 

14.80 

6.76 

6.89 

6.19 

12.37 

9.71 

11.58 

STD.DEVS. 

0.30 

0.23 

0.16 

0.12 

0.10 

0.11 

0.19 

0.13 

0.12 

X ERROR 

1.50 

1.44 

1.06 

1.72 

1.50 

1.72 

1.50 

1.32 

1.00 

GROUP 17 
AV ERAGES 

21 *33 

17.23 

15.50 

7.16 

7.28 

6.49 

13.08 

10.22 

12.11 

STO.DEVS. 

0*37 

0.28 

0.16 

0.13 

0. 14 

0.13 

0.24 

0.12 

0.11 

X ERROR 

1-75 

1.61 

1.02 

1.79 

1.94 

2.04 

1.81 

1.18 

0.89 

GROUP 18 
AVERAGES 

20*12 

16*38 

14.87 

6.72 

6.86, 

6.14 

12.42 

9.77 

11.63 

STO.DEVS, 

0*36 

0. 23 

0.13 

0.12 

0.14 

0.12 

0.20 

0.10 

0.11 

X ERROR 

1*80 

1.41 

0.88 

1.82 

2.06 

1.89 

1.56 

0.98 

0.93 

GROUP 19 
AVERAGES 

20*60 

16*65 

15.08 

6.90 

7.03 

6.28 

12.63 

9.91 

11.79 

STO.DEVS. 

0.39 

0*26 

0.22 

0.13 

0.15 

0.14 

0.20 

0.15 

0.17 

X ERRUR 

1*89 

1.59 

1.49 

1.90 

2.12 

2.15 

1.61 

1.51 

1.47 

GROUP 20 
AVERAGES 

20.06 

16.31 

14.77 

6.69 

6.84 

6.12 

12.37 

9.71 

11.55 

STD.DEVS. 

0*26 

0*20 

0.15 

0.10 

0.09 

0.08 

0.16 

0.10 

0.12 

X ERROR 

1.28 

1*23 

1.00 

1.45 

1.34 

1.35 

1.31 

1.07 

1.00 

GROUP 21 
AVERAGES 

20.54 

16.59 

15.00 

6.90 

7.01 

6.24 

12.60 

9.85 

11.73 

STD.DEVS. 

0.29 

0.25 

0.17 

0.10 

0.10 

0*11 

0.20 

0.14 

0.14 

S ERRUR 

1.43 

1.50 

1.15 

1.50 

1.49 

1.76 

1.58 

1.39 

1.17 

GROUP 22 
AVERAGES 

19.92 

16.41 

14.92 

6.57 

6.75 

6.13 

12.45 

9.82 

11.66 

STD.DEVS. 

0*11 

0.07 

0.07 

0.05 

0.05 

0.05 

0.05 

0.05 

0.06 

X ERROR 

0*54 

0*44 

0.47 

0.74 

0.75 

0.84 

0.40 

0.53 

0.54 

GROUP 23 
AVERAGES 

20*32 

16.47 

14.95 

6.78 

6.94 

6.21 

12.48 

9.81 

11.70 

STD.DEVS. 

0*26 

0. 19 

0.16 

0.08 

0.10 

0.09 

0.15 

0.12 

0.11 

X ERROR 

1.28 

1.18 

1.05 

1.13 

1.50 

1.53 

1.20 

1.22 

0.98 

GROUP 24 
AVERAGES 

19.50 

15.91 

14.54 

6.51 

6.65 

5.94 

12.07 

9.52 

11.38 

STD.DEVS. 

0.30 

0.21 

0.14 

0.13 

0.12 

0.09 

0.18 

0.09 

0. 12 

X ERROR 

1.52 

1.29 

0.95 

1.96 

1.73 

1.47 

1.46 

0.96 

1.04 




140. 

D2 {cont'd) 



17 

18 

19 

20 

21 

22 

VI 

V2 

V3 

GROUP 13 










AVERAGES 

20.29 

16.40 

14.79 

6.78 

6.92 

6.20 

12.45 

9.72 

11.57 

STO.DEVS. 

0.42 

0.30 

0.26 

0.17 

0.16 

0.11 

0.22 

0.18 

0.20 

X ERROR 

2.07 

1.80 

1.73 

2.44 

2.30 

1.81 

1.78 

1.87 

1.77 

GROUP 14 
AVERAGES 

21.29 

17.19 

15.52 

7.15 

7.26 

6.48 

13.05 

10.21 

12.13 

STD. DEVS. 

0.62 

0.43 

0.29 

0.26 

0.22 

0.21 

0.36 

0.22 

0.23 

X ERROR 

2.92 

2.53 

1.84 

3.68 

3.06 

3.19 

2.75 

2.15 

1 . 86 

GROUP 15 
AVERAGES 

20.48 

16.59 

14.96 

6. 86 

6.98 

6.24 

12.60 

9.84 

11.69 

STD. DEVS. 

0.23 

0.18 

0.20 

0.10 

0.10 

0.06 

0. 14 

0.12 

0.18 

X ERROR 

1.13 

l.ll 

1.31 

1.41 

1.37 

1.04 

1.14 

1.21 

1.51 

GROUP 16 
AVERAGES 

20.22 

16.32 

14.80 

6.76 

6.89 

6.19 

12.37 

9.71 

11.58 

STD. DEVS. 

0.30 

0.23 

0.16 

0.12 

0.10 

0.11 

0.19 

0.13 

0.12 

X ERROR 

1.50 

1.44 

1.06 

1.72 

1.50 

1.72 

1.50 

1.32 

1.00 

GROUP 17 
AVERAGES 

21.33 

17.23 

15.50 

7.16 

7.28 

6.49 

13.08 

10.22 

12.11 

STD. DEVS. 

0.37 

0.28 

0.16 

0.13 

0.14 

0.13 

0.24 

0.12 

0.11 

X ERROR 

1.75 

1.61 

1.02 

1.79 

1.94 

2.04 

1.81 

1.18 

0.89 

GROUP 18 
AVERAGES 

20.12 

16.38 

14.87 

6.72 

6.86 

6.14 

12.42 

9.77 

11.63 

STO.DEVS. 

0.36 

0.23 

0.13 

0.12 

0.14 

0.12 

0.20 

0.10 

0.11 

X ERROR 

1.80 

1.41 

0.88 

1.82 

2.06 

1.89 

1.58 

0.98 

0.93 

GROUP 19 
AVERAGES 

20.60 

16.65 

15.08 

6.90 

7.03 

6.28 

12.63 

9.91 

11.79 

STD.DEVS. 

0.39 

0.26 

0.22 

0.13 

0.15 

0.14 

0.20 

0.15 

0.17 

X ERROR 

1.89 

1.59 

1.49 

1.90 

2.12 

2.15 

1.61 

1.51 

1.47 

GROUP 20 
AVERAGES 

20.06 

16.31 

14.77 

6.69 

6.84 

6.12 

12.37 

9.71 

11.55 

STD.DEVS. 

0.26 

0.20 

0.15 

0.10 

0.09 

0.08 

0.16 

0.10 

0.12 

X ERROR 

1.28 

1.23 

1.00 

1.45 

1.34 

1.35 

1.31 

1.07 

1.00 

GROUP 21 
AVERAGES 

20.54 

16.59 

15.00 

6.90 

7.01 

6.24 

12.60 

9.85 

11.73 

STD.DEVS. 

0.29 

0.25 

0.17 

0.10 

0.10 

0.11 

0.20 

0.14 

0.14 

X ERROR 

1-43 

1.50 

1.15 

1.50 

1.49 

1.76 

1*58 

1.39 

1.17 

GROUP 22 
AVERAGES 

19.92 

16.41 

14.92 

6.57 

6.75 

6.13 

12.45 

9.82 

11.66 

STD.DEVS. 

0.11 

0.07 

0.07 

0.05 

0.05 

0.05 

0.05 

0.05 

0.06 

X ERROR 

0.54 

0.44 

0.47 

0.74 

0.75 

0.84 

0.40 

0.53 

0.54 

GROUP 23 
AVERAGES 

20.32 

16.47 

14.95 

6.78 

6.94 

6.21 

12.48 

9.81 

11.70 

STD.DEVS. 

0.26 

0.19 

0.16 

0.08 

0.10 

0.09 

0.15 

0.12 

0.11 

X ERROR 

1.28 

1.18 

1.05 

1.13 

1.50 

1.53 

1.20 

1.22 

0.98 

GROUP 24 
AVERAGES 

19.50 

15.91 

14.54 

6.51 

6.65 

5.94 

12.07 

9.52 

11.38 

STD.DEVS. 

0.30 

0.21 

0.14 

0.13 

0.12 

0.09 

0.18 

0.09 

0. 12 

X ERROR 

1.52 

1.29 

0.95 

1.96 

1.73 

1.47 

1.46 

0.96 

1.04 



140 . 

B2 {qont'd) 




17 

18 

19 

20 

21 

22 

VI 

V2 

V3 

GROUP 13 








~ L ' 


AVERAGES 

20.29 

16.40 

14.79 

6.78 

6.92 

6.20 

12.45 

9.72 

11.57 

STD.DEVS. 

0.42 

0.30 

0.26 

0.17 

0.16 

0.11 

0.22 

0.18 

0.20 

X ERROR 

2.07 

1.80 

1.73 

2.44 

2.30 

1.81 

1.78 

1.87 

1.77 

GROUP 14 
AVERAGES 

21.29 

17.19 

15.52 

7.15 

1 7.26 

6.48 

13.05 

10.21 

12.13 

STD.DEVS. 

0.62 

0.43 

0.29 

0.26 

0.22 

0.21 

0.36 

0.22 

0.23 

* ERROR 

2.92 

2.53 

1.84 

3.68 

3.06 

3.19 

2.75 

2.15 

1.86 

GROUP IS 
AVERAGES 

20.48 

16.59 

14.96 

6.86 

6.98 

6.24 

12.60 

9.84 

1 1.69 

STD.OEVS. 

0.23 

0.18 

0.20 

0.10 

0.10 

0.06 

0.14 

0.12 

0.18 

X ERROR 

1.13 

1.11 

1.31 

1.41 

1.37 

1.04 

1.14 

1.21 

1.51 

GROUP 16 
AVERAGES 

20.22 

16.32 

14.80 

6.76 

6.89 

6.19 

12.37 

9.71 

11.58 

STD.DEVS. 

0.30 

0.23 

0.16 

0.12 

0.10 

0.11 

0.19 

0.13 

0.12 

X ERROR 

1.50 

1.44 

1.06 

1.72 

1.50 

1.72 

1.50 

1.32 

1.00 

GROUP 17 
AV ERAGES 

21.33 

17.23 

15.50 

7.16 

7.28 

6.49 

13.08 

10.22 

12.11 

STD.OEVS. 

0.37 

0.28 

0.16 

0.13 

0.14 

0.13 

0.24 

0.12 

0.11 

X ERROR 

1.75 

1.61 

1.02 

1.79 

1.94 

2.04 

1.81 

1.18 

0.89 

GROUP 18 
AVERAGES 

20.12 

16.38 

14.87 

6.72 

6.86. 

6.14 

12.42 

9.77 

11.63 

STD.DEVS. 

0.36 

0. 23 

0.13 

0.12 

0.14 

0.12 

0.20 

0.10 

0.11 

X ERROR 

1.80 

1.41 

0.88 

1.82 

2.06 

1.89 

1.58 

0.93 

0.93 

GROUP 19 
AVERAGES 

20.60 

16.65 

15.08 

6.90 

7.03 

6.28 

12.63 

9.91 

11.79 

STD.DEVS. 

0.39 

0.26 

0.22 

0.13 

0.15 

0.14 

0.20 

0.15 

0.17 

X ERROR 

1.89 

1.59 

1.49 

1.90 

2.12 

2.15 

1.61 

1.51 

1.47 

GROUP 20 
AVERAGES 

20.06 

16.31 

14.77 

6.69 

6.84 

6.12 

12.37 

9.71 

11.55 

STD.DEVS. 

0.26 

0.20 

0.15 

0.10 

0.09 

0.08 

0.16 

0.10 

0.12 

X ERROR 

1.28 

1.23 

1.00 

1.45 

1.34 

1.35 

1.31 

1.07 

1.00 

GROUP 21 
AVERAGES 

20.54 

16.59 

15.00 

6.90 

7.01 

6.24 

12.60 

9.85 

11.73 

STD.OEVS. 

0.29 

0.25 

0.17 

0.10 

0.10 

0.11 

0.20 

0.14 

0.14 

X ERROR 

1.43 

1.50 

1.15 

1.50 

1.49 

1.76 

1.58 

1.39 

1.17 

GROUP 22 
AVERAGES 

19.92 

16.41 

14.92 

6.57 

6.75 

6.13 

12.45 

9.82 

11.66 

STD.OEVS. 

0.11 

0.07 

0.0? 

0.05 

0.05 

0.05 

0.05 

0.05 

0.06 

X ERROR 

0.54 

0.44 

0.47 

0.74 

0.75 

0.84 

0.40 

0.53 

0.54 

GROUP 23 
AVERAGES 

20.32 

16.47 

14.95 

6.78 

6.94 

6.21 

12.48 

9.81 

11.70 

STO.DEVS. 

0.26 

0.19 

0.16 

0.08 

0.10 

0.09 

0.15 

0.12 

0.11 

X ERROR 

1.28 

1.18 

1.05 

1.13 

1.50 

1.53 

1.20 

1.22 

0.98 

GROUP 24 
AVERAGES 

19.50 

15.91 

14.54 

6.51 

6.65 

5.94 

12.07 

9.52 

11.38 

STD. 3 tVS. 

0.30 

0.21 

0.14 

0.13 

0.12 

0.09 

0.18 

0.09 

0. 12 

X ERROR 

1.52 

1.29 

0.95 

1.96 

1.73 

1.47 

1.46 

0.96 

1.04 
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141. 
D2 (cont'd) 



17 

18 

19 

20 

21 

22 

VI 

V2 

V3 

GROUP 25 










AVERAGES 

20.38 

16.53 

14.96 

6.81 

6.96 

6.20 

12.54 

9.83 

11.70 

STO.OEVS. 

0.48 

0.35 

0.22 

0. 17 

0.18 

0.14 

0.29 

0. 19 

0.17 

* ERROR 

2.37 

2.14 

1.48 

2.51 

2.66 

2.2 5 

2.29 

1.96 

1.43 

GROUP 26 

AVERAGES 

20*41 

16.55 

15.01 

6.80 

6.97 

6.24 

12.55 

9.85 

11.74 

STO.OEVS. 

0*22 

0.15 

0.13 

0.09 

0.09 

0.06 

0.13 

0.08 

0.12 

% ERROR 
GROUP 27 

1.06 

0.92 

0.67 

1.37 

1.25 

0.98 

1.01 

0.83 

0.98 

AVERAGES 

19.76 

16.11 

14.67 

6.58 

6.74 

6.03 

12.21 

9.63 

11.48 

STO.OEVS. 

0.19 

0.15 

0.09 

0.09 

0.08 

0.13 

0.12 

0.06 

0.07 

% ERROR 

0*98 

0.91 

0.58 

1*36 

1.13 

2.12 

0.96 

0.67 

0.59 

GROUP 28 

AVERAGES 

19.95 

16*28 

14.82 

6.70 

6.80 

6*04 

12.35 

9.74 

11.59 

STO.OEVS. 

0*22 

0.14 

0.10 

0.07 

0.10 

0*14 

0*11 

0.06 

0.09 

X ERROR 
GROUP 29 

1*12 

0*86 

0.69 

1.00 

1.41 

2*31 

0*86 

0.67 

0.74 

AVERAGES 

20.04 

16.34 

14.80 

6.65 

6.83 

6*14 

12*40 

9.73 

11.57 

STO.OEVS. 

0.22 

0*19 

0.18 

0.08 

0.08 

0*09 

0.15 

0.10 

0.16 

X ERROR 

1.08 

1.13 

1.22 

1.27 

1.18 

1*41 

1*24 

1.07 

1.39 

GROUP 30 

AVERAGES 

20.64 

16.94 

15.36 

6. 86 

7.00 

6*33 

12.85 

10.12 

12.00 

STD. DEVS. 

0.08 

0.06 

0.05 

0.05 

0.04 

0*04 

0*06 

0.05 

0.06 

X ERROR 
GROUP 31 

0*37 

0.35 

0.35 

0.66 

0.57 

0*66 

0.46 

0.50 

0.48 

AVERAGES 

20.65 

16.93 

15.34 

6.86 

7.00 

6.33 

12*84 

10.11 

11.98 

STO.OEVS. 

0*13 

0.09 

0.09 

0.05 

0.06 

0.04 

0*08 

0.06 

0.08 

S ERROR 
GROUP AO 

0.62 

0*52 

0.58 

0*79 

0.81 

0.67 

0*62 

0.63 

0.65 

AVERAGES 

21.33 

17.23 

15.48 

7.13 

7.29 

6.51 

13*07 

10.21 

12.09 

STO.OEVS. 

0*43 

0* 32 

0.31 

0.16 

0.16 

0.12 

0.25 

0.20 

0.25 

X ERROR 

2*02 

1.87 

1.98 

2.24 

2.15 

1.88 

1*93 

2.00 

2.05 

GROUP 41 

AVERAGES 

21.22 

17. 17 

15.46 

7.08 

7.25 

6.49 

13*03 

10.18 

12.08 

STO.OEVS. 

0*48 

0.33 

0.29 

0.19 

0.17 

0*13 

0*26 

0. 19 

0.23 

X ERROR 
GROUP 42 

2.24 

1.94 

1.87 

2.70 

2.38 

2*06 

1.98 

1.88 

1.90 

AVERAGES 

20*49 

16.57 

14.99 

6.87 

7.01 

6*25 

12.57 

9.85 

11.73 

STO.OEVS. 

0*51 

0.40 

0.34 

0.19 

0.17 

0*16 

0*30 

0.24 

0.26 

X ERROR 

2.50 

2.41 

2.29 

2.73 

2.48 

2*57 

2.42 

2.43 

2.25 

GROUP 43 

AV ERAGES 

22.06 

17.70 

15.82 

7.40 

7.53 

6.74 

13.44 

10.43 

12.36 

STO.OEVS. 

0.66 

0.53 

0.40 

0.24 

0.22 

0.21 

0.42 

0.28 

0.31 

f ERROR 

2.99 

2.99 

2.53 

3.24 

2.89 

3.13 

3.15 

2.67 

2.49 

GROUP 44 

AVERAGES 

20.84 

16.77 

15.06 

6*96 

7*13 

6*39 

12.72 

9.93 

11.76 

STO.OEVS. 

0.33 

0.27 

0.20 

0.11 

0*08 

0.14 

0.20 

0. 16 

0.14 

X ERROR 

1.56 

1.61 

1.33 

1.59 

1*18 

2.14 

1.56 

1.64 

1.17 
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17 

18 

19 

20 

21 

22 

VI 

. V2 

V3 

GROUP 45 
AVERAGES 

20.51 

16.53 

15. L2 

6.83 

7.01 

6.29 

12.51 

9.87 

11.86 

STD.OEVS. 

o.sa 

0.39 

0.46 

0.30 

0.23 

0.17 

0.32 

0.25 

0.38 

X ERROR 

2.83 

2.38 

3.03 

4.39 

3.31 

2.71 

2.56 

2.50 

3.25 

GROUP 46 
AVERAGES 

21.38 

17.21 

15.52 

7.12 

7.31 

6.56 

13.05 

10.18 

12.14 

STD.OEVS. 

0.64 

0.45 

0.39 

0.31 

0.23 

0.16 

0.35 

0.28 

0.32 

1 ERROR 

3.00 

2.63 

2.52 

4.35 

3.09 

2.47 

2.69 

2.77 

2.61 

GROUP 50 
AVERAGES 

21.28 

17.19 

15.18 

7.09 

7.10 

6.68 

13.03 

9.99 

11.81 

STD.OEVS. 

0.08 

0.14 

0.22 

0. 03 

0.03 

0.05 

0.09 

0.14 

O.ZO 

X ERROR 

0.38 

0.79 

1.47 

0.39 

0.48 

0.75 

0.72 

1.37 

1.69 

GROUP 5 i 
AVERAGES 

21.66 

17.41 

15.39 

7.21 

7.35 

6.77 

13.21 

10.17 

11.95 

STD. DEVS. 

0.22 

0.18 

0.59 

0.10 

0.08 

0.03 

0.16 

0.07 

0.64 

X ERROR 

1.01 

1.01 

3.86 

1.37 

1.06 

0.42 

1.23 

0.70 

5.38 

GROUP 52 
AVERAGES 

21.48 

17.49 

15.75 

7.06 

7.24 

6.78 

13.28 

10.24 

12.32 

STD. DEVS. 

0.02 

0.01 

0.02 

0.03 

0.02 

0.03 

0.03 

0.04 

0.04 

% ERROR 

0.10 

0.06 

0.12 

0.44 

0.32 

0.49 

0.24 

0.40 

0.31 

GROUP 53 
AVERAGES 

17.50 

14.38 

L3.19 

5.75 

5.96 

5.41 

10.89 

8.62 

10.33 

STD.OEVS. 

0.05 

0.03 

0.02 

0.03 

0.03 

0.03 

0.04 

0.04 

0.05 

X ERROR 

0.30 

0.23 

0.15 

0.51 

0.52 

0.52 

0.39 

0.46 

0.49 

GROUP 54 
AVERAGES 

20.22 

16.44 

14.89 

6.75 

6.89 

6.18 

12.47 

9.78 

11.64 

STO.DEVS. 

0.37 

0.28 

0.22 

0.14 

0.14 

0.12 

0.22 

0.16 

0.17 

X ERROR 

1.84 

1.68 

1.49 

2.03 

2.06 

2.00 

1.76 

1.69 

1.50 
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O < 4 *# 

3 « 4 •# 

3«* 4 *# 

3 < 4 ** 


GROUP ll GROUP 23 GROUP O 

AVERAGES 1.2346 1.3671 1.1074 AVERAGES 1.2341 1.3591 1.1013 AVERAGES 1.2460 1.3941 1.1168 
S1U.0EVS. 0.0065 0.0270 0.0182 SIO.DEVS. 0.0056 0.0085 0.0040 SIU.OEVS. 0.0100 0.0214 0.0142 
X ERROR 0.52 L.97 1.65 % ERROR 0.46 0.62 0.37 t ERROR 0.61 1.53 1.27 
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COMPARISON OF AIRBORNE INFRARED SPECl'RAL EMITTANCE AND 
RADAR SCAT1ERCMLTLR DATA FROM PISGAH CRATER LAVA FLOWS 


R. J. P. Lyon 

Stanford University 
Stanford, California 

SUMMARY 


The olivine basalt lava fields of Pisgah Crater, 35 miles KSE of Barstow, California, are one 
of the very few areas which have boon studied by more than one sensor. In fact, a big problem in 
evaluating vai*ious remote sensing sy steins is that so rarely have they been viewing the same targets, 
let alone .from the same altitude aim 1 uok angles. 

In NASA/MSC Mission 10 8 the infrared spectrornetcr/radiometer instrumentation was flown twice 
down the same trsree-miJe long flight line, as was used for the 2 cm-band radar scatterometer in 
Mission 21, anu imported at the Fifth Symposium on Remote Sensing of Environment. These two, 
non- imaging systems moasui'e uiiTer*ent parameters - spectral emittance for the infrared, and gonio- 
metric radar backscatier for tne scatterorr.etur. In addition, cross-track width of the "ground patch" 
(or footprint), of the Ik units is 7 iMllirauians and the radar is ^ milliradians . Despite these 
obvious differences, cat a from both units can be used to. arrive at similar classifications for the 
geological materials at this site - lava flows of three types, lava cinders of Pisgah Crater, dry 
sediments of La vie Lake, and several types of Older and Younger alluviums in the desert outwash fans 
surrounding the crater arid flows. 


1449 


APPENDIX TABLE D5 


SEQUENTIAL SUCCESS BY ROCK TYPE FOR EACH METHOD USING BMDQ7M 


A. MIXED TERRAIN (D, A, P) 


B. ALL LAVAS (I, II, III) 


1. Spectral Emittance 
(50 central X's) 


1. Spectral Emittance 
(50 central X's) 




D 

A 

p 


I 

II 

III 

Step 

1 

70 

93 

35 

Step 1 

66 

27 

54 


2 

68 

100 

48 

2 

78 

35 

51 


3 

68 

96 

52 

3 

69 

48 

51 


4 

81 

100 

77 

4 

66 

55 

49 


6 

81 

100 

78 

6 

72 

55 

52 


7 

84 

100 

81 

7 

66 

57 

55 


9 

85 

100 

87 

9 

69 

58 

57 


2. 

System 

Response 

2. 

System Response 



a. Vincent's 

- 3 


a. Vincent' 

s - 3 

Step 

1 

87 

100 

87 

Step 1 

52 

61 

63 


2 

87 

100 

97 

2 

61 

57 

69 


3 

96 

100 

97 

3 

65 

55 

66 



b. MSDS - 6 



b. MSDS 

- 6 


Step 

1 

89 

100 

84 

Step 1 

48 

26 

67 


2 

100 

100 

84 

2 

65 

52 

69 


3 

100 

100 

97 

3 

81 

55 

66 


4 

100 

100 

97 

4 

75 

47 

65 


6 

100 

100 

97 

6 

74 

59 

63 


3. Ratios 3. Ratios 


a. Vincent's - 3 a. Vincent's - 3 


Step 1 

87 

100 

71 

Step 1 

58 

11 

63 

2 

83 

100 

71 

2 

56 

44 

57 

b. 

MSDS 

- 15 


b. 

MSDS 

- 15 


Step 1 

51 

100 

32 

Step 1 

22 

15 

62 

2 

72 

100 

68 

2 

68 

26 

58 


149 . 



